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Probing the collective spin state of an ensemble of atoms may provide a means to reduce heating 
via the photon recoil associated with the measurement and provide a robust, scalable route for 
preparing highly entangled states with spectroscopic sensitivity below the standard quantum limit 
for coherent spin states. The collective probing relies on obtaining a very large optical depth that 
can be effectively increased by placing the ensemble within an optical cavity such that the probe light 
passes many times through the ensemble. Here we provide fundamental expressions for such cavity- 
aided non-demolition measurements. We also provide experimental examples related to previous 
work [I] using quantum non-demolition measurements to prepare conditionally spin-squeezed states 
in a large ensemble of nearly 10 6 87 Rb atoms. Fundamental limits on spectroscopic enhancements 
in 87 Rb are considered. 

PACS numbers: 42.50.Dv, 42.50.Pq, 06.20.-f, 42.50.Lc 



I. INTRODUCTION 

High-resolution measurements of the populations of 
two-level systems are key for realizing high precision 
atomic sensors such as atomic clocks, magnetometers, 
and atom-based electric field, rotation, and inertial sen- 
sors [2]- Further, developing non-demolition, high- 
resolution measurement techniques to create and/or de- 
tect entangled states is a promising route to enhanced 
sensor precision or bandwidth [TJ l3TfnT] . 

In trapped neutral atom ensembles, non-demolition 
measurements that do not cause atom loss from the trap 
could lead to significant advances in the repetition rates 
of sensors [13], allowing them to operate closer to the 
regime of ion-based sensors in which the ions can be 
stored over many repeated measurement cycles [2] . Fur- 
thermore, a quantum non-demolition measurement that 
also preserves quantum coherence can prepare condition- 
ally spin-squeezed states with spectroscopic sensitivity 
below the Standard Quantum Limit (SQL) A#sql = 
l/V^V that arises from the quantum projection noise of 
N independent atoms [T5] . 

The results described in this work are relevant to the 
recently demonstrated approaches for generating entan- 
gled states in large ensembles using many diverse ap- 
proaches, including quantum non-demolition measure- 
ments [TJiniE], one-axis twisting schemes based on cavity- 
mediated interactions [16] . and direct collisional interac- 
tions that generate one-axis squeezing [5j |9] or para- 
metric pair generation [H]. In all of these cases, a 
low noise readout such as the approach described in this 
paper is always required to actually exploit the enhanced 
phase-sensing properties of these states. 

In particular, there has been substantial recent interest 
in developing non-demolition readout schemes for large 
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FIG. 1. (Color online) (a) Relevant energy levels. The 
pseudo-spin 1/2 system typically comprises two metastable 
states |4-), which are utilized in atomic sensors and clocks. 
The number of atoms Nf in |t) modifies the cavity resonance 
frequency. Initially, we will assume the optically excited state 
|e) does not couple to \\) because of dipole selection rules or 
because the coupling is highly non-resonant S c <C ujm, where 
ujm is the frequency separation between | f ) and | J,) . (b) 
Transmitted and/or reflected probe light is monitored to de- 
termine the cavity resonance frequencies and hence the num- 
ber of atoms in |f) . An applied NMR-like microwave rotation 
can swap the populations between ground states so that a 
subsequent measurement can also determine the |^) popula- 
tion N± . The ensemble of atoms are trapped in an intracavity 
optical lattice (blue-green). The total cavity power decay rate 
is k. Single atom spontaneous decay from |e) at rate T leads 
to photon recoil heating and single-atom wavefunction col- 
lapse. The goal then is to extract collective information from 
the probe mode more rapidly than the undesired single-atom 
photon scattering into free space. 



laser-cooled and quantum-degenerate neutral atomic en- 
sembles consisting of roug hly 10 3 to 10 7 atoms [TJ 
HTl [TBI [T8H2fJ] . It is well known that significant improve- 
ments in readout sensitivity can be achieved by optically 
probing ensembles in free space along directions of high 
resonant optical depth. This approach has been exten- 
sively analyzed theoretically , and experimentally 
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FIG. 2. (Color online) Visualization of a non-demolition (ND) 
(upper right) and a quantum non-demolition (QND) (lower 
right) measurement. A coherent spin state (CSS) for an en- 
semble of N spin- 1/2 atoms prior to measurement is repre- 
sented as a Bloch vector (red arrow) of length N/2. The 
quantum noise in the orientation of the Bloch vector is visu- 
alized as a quasi-probability distribution (red/yellow region) 
perpendicular to the Bloch vector, with rms opening angle 
at the Standard Quantum Limit (SQL) A#sql = 1/yN. A 
non-demolition measurement in this paper will refer to a mea- 
surement of the state's z spin projection J z = (Nf — N^)/2 
with an rms imprecision AJ Z < AJ Zj css = V^V /2, and with 
the majority of the atoms remaining trapped after the mea- 
surement. For example, in the ND measurement visualized on 
the right, the measurement imprecision is AJ Z = 0, but after 
the measurement the atoms are described by a product state 
of atoms in spin up and down due to single-atom state infor- 
mation gained by the environment via free space scattering of 
light. We define a quantum non-demolition measurement via 
the additional requirement that a sufficiently large number of 
atoms remain in a coherent superposition of spin up and down 
such that the resulting state, conditioned on the measure- 
ment outcome, has a polar angle uncertainty A9 < A#sql- 
The definition of quantum non-demolition employed here is 
related to, but less restrictive than that of Ref. [17] . The 
lower right QND example visualizes the conditional state as 
a product state of a squeezed state and the atoms that have 
been collapsed into spin up and down. 
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More recently, the technique of free space probing of 
high optical depth samples has been extended to using 
optical cavities to effectively increase the optical depth of 
the atomic ensemble by passing the probe light through 
the same sample many times as shown in Fig. [j] [Tl 171 HT1- 
146] . In most of these experiments and proposals, the 
cavity is far-detuned from the optical transition that was 
probed. Probing in resonant regime [T] is an exception 
rather than the norm. In principle, the cavity detuning 
8 C can be chosen almost arbitrarily. Therefore, a natural 
question to ask is: How does cavity detuning affect both 
the fundamental and technical atomic population mea- 
surement resolution or the degree of spin-squeezing for a 



given cavity geometry and cavity finesse? 

To answer the question posed above, we provide de- 
tailed expressions for the fundamental scalings for prob- 
ing an atomic ensemble using an optical cavity that 
smoothly connects the resonant to the far-detuned prob- 
ing regime. We apply our results to first estimate the 
amount of photon-recoil heating of the ensemble when 
the cavity-aided measurement has an imprecision at the 
quantum projection noise level. The average number of 
photon recoils per atom sets the degree to which the mea- 
surement can be considered non-demolition. This analy- 
sis is then extended to estimate bounds on the degree of 
conditional spin squeezing (see Fig.[2]for details) that can 
be obtained using 87 Rb atoms in a cavity (see Fig. 
[H [7J. We show that the fundamental limitations are 
set by the collective cooperatively parameter NC, the 
probe detection quantum efficiency g,and atomic prop- 
erties alone. The collective cooperativity parameter NC 
plays a similar role to the resonant optical depth of atoms 
in free space, where N is the number of atoms in the 
probe volume and C is the single particle cooperativ- 
ity parameter [37]. To illustrate these concepts we also 
consider technical advantages and experimental methods 
used in Ref. [T] for the special case that the cavity is 
tuned to resonance with the optically excited state. 

This paper is organized as follows. In Sec. [TTJ we begin 
with a review of the properties of the coupled atoms- 
cavity system including dissipation. This review also 
provides definitions and notation used throughout the 
paper. 

In Sec. |III[ we derive the quantum-limited signal-to- 
noise ratio for non-demolition measurements of atomic 
populations considered as a function of the cavity detun- 
ing from atomic resonance (S c in Fig. [l^.). We identify 
three different probing regimes, a resonant regime and 
two detuned regimes separated by a critical detuning 5° . 



In Sec. IV we focus on our experimental system [T], in 
which the atoms are probed using a cavity resonant with 
the atomic optically excited state 5 C = 0. We include dis- 
cussions of technical advantages and limitations associ- 
ated with operating in this regime. Experimental details 
are provided for two probing schemes: a Pound-Drever- 
Hall type measurement and a phase-coherent demodula- 
tion of both the reflected and transmitted cavity fields 
into two orthogonal quadratures. We show a high degree 
of noise cancellation and experimentally demonstrate ac 
Stark shift cancellation enabled by simultaneous probing 
of the two resonances of the vacuum Rabi splitting. 

In Sec. [Vj we address quantum back-action effects due 
to probe-induced spin flips on estimates of atomic popu- 
lations in a simple three-level model. We show how the 
optimal measurement resolution can be achieved by bal- 
ancing between noise added by spin flips and averaging 
down the probe's vacuum or photon shot noise. 

In Sec. VI we consider the limits set on coherence 
preservation. Coherence is lost due to wavefunction col- 
lapse into spin up or down driven by the same probe 
induced free space scattering that also causes photon re- 
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coil heating. We then obtain the optimal spectroscopic 
enhancement, the figure of merit for a quantum non- 
demolition measurement, as a function of the spin flip 
probability p. 

In Sec. |VII[ we apply the results of Sec. [Vl] to two 
concrete examples in 87 Rb: generating conditional spin 
squeezing first using a non-cycling optical transition, and 
then a cycling optical transition. Here, we demonstrate 
the key role of the ratio of the ground state hyperfine 
splitting Whf to the optical transition width T for deter- 
mining scalings and fundamental limits on conditional 
spin squeezing. 

Finally, in Sec. VIII[ we theoretically show that fol- 
lowing a quantum non-demolition (QND) measurement 
near the projection noise level, a coherent re-initialization 
scheme can be used to reduce photon recoil heating as- 
sociated with re-preparing the atomic ensemble for the 
next experimental cycle. 



II. COUPLED ATOM CAVITY MODES 

To begin, we provide a brief review of the open coupled 
atoms-cavity system with the goal of both providing a 
framework for understanding the experimental work and 
for explicitly enumerating the assumptions made to re- 
duce this system to a classical two- mode system [3BH5U] . 
The dynamics of the system under a classical drive and 
dissipation are then studied with the goal of obtaining 
the full complex response of the reflected and transmit- 
ted cavity field. Finally, a discussion of the probe signal- 
to-noise sets the stage for addressing measurement reso- 
lution at the projection noise level in Sec. |III| 
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A. Linearized response 

We consider atoms with two ground states |f) and |T) 
whose populations we wish to estimate [see Fig. [TJb)]. 
Atoms in |f ) can interact with a cavity mode by absorb- 
ing a cavity photon and being promoted to an optically 
excited state |e). On the other hand, atoms in ||) are 
assumed to not interact with the cavity mode because 
of dipole selection rules, a large energy splitting between 
the ground states, or otherwise. We note that a coher- 
ence can exist between |f) and as would be required 
for sensing a quantum phase, but this coherence is not 
important for the discussion in this section. In the ab- 
sence of dissipation, the coupled atoms-cavity system is 
well described by the Tavis-Cummings Hamiltonian: 
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H = h5 c c) c + h gi 



a i+ c 
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(1) 



The Hamiltonian is written in a frame rotating at the 
atomic transition I'f) — > |e) angular frequency ujf e . The 
parameter g t sets the coupling strength between the cav- 
ity field and the ith atom. There are N atoms total. 
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FIG. 3. (Color online) (a) Graphical representation of the lin- 
earization achieved via the Holstein-PrimakofF approximation. 
The atomic sub-system \f) and |e) is described by a Bloch vec- 
tor of length (Nf + JV c )/2. In the regime in which the probe 
light only weakly excites the atoms such that N e <C N^, the 
two-level system can be described by its projection onto a 2D 
plane equivalent to that describing a light field, (b) With this 
approximation, the collective atomic mode may be treated as 
an equivalent cavity mode (lower purple mode) whose cou- 
pling to the actual cavity (upper red mode) is governed by 
a partially reflecting mirror (center) described by a field cou- 
pling rate constant (y/Nfg) that depends on the number of 
atoms in state jf) • The physical mirrors have transmission co- 
efficients described by the field coupling rates ^/k^ and s/k2, 
and internal cavity losses are described by %/kl, such that 
the total power decay rate is k — Ki + k.2 + Kl- Decay of the 
atoms by emission of a photon into free space is described by 
the field transmission coefficient VT that is the sum of the 
field scattering into free space modes by the atoms in f) . 
Because the scattered modes are distinguishable (expanded 
view of free space scattering port), it is possible to tell which 
atoms are in | f) from the free-space-scattered photons, de- 
stroying any coherent superposition between |f) and another 
state |4_) (not shown) that may have been prepared for sensing 
a quantum phase. 



The operators & and c are the creation and annihilation 
operators for the cavity field with commutation relations 
[c, cH = 1 and photon number operator M c = &c. The 
cavity detuning S c = cu c — ujf e is the difference between 
the empty cavity frequency oj c and the atomic transition 
frequency. Finally, the raising and lowering operators for 
the ith atom are &i + — |ej)(ti|, and <7j_ = |ti)(ej| re- 
spectively. The number of atoms in each state is given 
by the collective projection operators Nf = J2iLi lti)(t»|j 



N i = Eili \U){U\, an d N e = J2"=i |e l )(e i |. For brevity, 
we use the following abbreviations throughout this paper: 



A? 



t = \ N t)> N i = yn / 

Although the atoms may in general exist in a superpo- 
sition of spin up and spin down, in the following analysis 
we will consider the atoms to be in a definite eigenstate 
Nf of the Nf operator. We will then reintroduce the fluc- 
tuations in the operator Nf for atoms in a superposition 



N\ ) and N P 
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of |t) and |4_) in terms of the rms projection noise about 



the mean value AJV- 



t 



(JV t - (N t ))' 



We assume that the system is driven weakly so that 
the mean number of atoms in the optically excited state 
|e) is a small fraction of the total number of atoms in 
|t), such that N e /N^ <C 1. We also assume that the 
coupling is uniform gi = g (with appropriate renormal- 
ization of g and N in the case of inhomogeneous cou- 
pling provided in Ref. pQ). In the weak excitation limit, 
one can make the Holstein-Primakoff approximation [51] 
that essentially replaces the three-dimensional Bloch vec- 
tor describing the two-level system composed of states |e) 
and |t) by its vector projection onto the two dimensional 
x — y plane [Fig. |3ja)]. When no probe light is applied, 
the system relaxes into its ground state orientation — z. 
When the system is weakly excited by probe light, and a 
small fraction of the atoms is promoted to the optically 
excited state, the lowest order effect on the Bloch vec- 
tor is a small angular deflection of the Bloch vector away 
from — z. 

In this approximation, one replaces the atomic ladder 
operators with effective creation and annihilation oper- 

3 y/NftfJl 



ators as Y^i ® 



a 



and Ef &i 



N e « 
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In the limit of low atomic saturation 



duce to 



N+, the creation and annhilation operators re- 



X)f £»+> and 



a 



Ei °"i 



The 



number of atoms in |e) is given by the number operator 
N e = a f a. 

The resulting Hamiltonian describes the coherent cou- 
pling of two ring cavities through a beam-splitter whose 
transmission coefficient increases with the number of 
atoms in |t), i.e., 



H = hSrc'c ■ 



l^g (ac f + a)c) 



(2) 



This model is shown in Fig. |3jb). The coupled Heisen- 
berg equations of motion for the cavity and atomic oper- 
ators are 



dt 
d(a) 
dt 



-i8 c (c) - lyfN^g (a) , 



T3 ( 



(3) 



We can solve for the eigenfrequencies of the coupled 
equations in terms of the cavity detuning S c and the col- 
lective vacuum Rabi frequency ft^- = yjN^2g using 



5r.±JSi 



(4) 



The collective vacuum Rabi splitting O-j- sets the differ- 
ence in the normal mode frequencies lj + — w_ at zero 
detuning 5 C — 0. The normal mode frequencies in the 
laboratory frame are io^ e + lo±. A useful relation for 
simplifying expressions involving the mode frequencies is 
cj + w_ = -0,1/4. 
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FIG. 4. (Color online) The normal modes of the coupled 
atoms-cavity system h± exhibit an avoided crossing in the 
transmission spectrum at 8 C = 0, where the size of the collec- 
tive vacuum Rabi splitting is Q,^ = \fW^2g. 



Any arbitrary superposition of atomic and cavity am- 
plitudes can be written in terms of complex normal 
mode amplitudes B± and unit vectors h± as ((c) ,(&)) = 



B + b + + e _i "-*S_b_. At resonance S c = 0, the 



normal mode vectors are h± = (l,±l)/\/2, and the nor- 
mal modes beat against one another at relative frequency 
Q^. For example, if at t = there are no atoms in the 
excited state, then B± = l/v2- The light in the cavity 
will Rabi flop (or be converted) into purely excited-state 
atoms after a time 7r/fi^, and then Rabi flop back into 
light in the cavity after 2ir/Vl^. 

More generally, the unit vectors can be expressed in 
terms of a mode-mixing angle 9 rn = arctan(w_/w + ) 
as b + = (cos# m , sin0 m ), and b_ = (— sin0 m , cos m ). 
When S c S> fi-t-, the b + normal mode is predominantly 
cavity-like, while the b_ normal mode is predominantly 
atom-like. The converse is true when 5 C <C —flf. This 
normal mode structure is demonstrated in Fig. [4] with 
an experimental measurment of the transmitted probe 
power through the cavity as a function of the probe laser 
frequency and the cavity detuning S c . 



B. Driven and damped dynamics 

The Heisenberg equations of motion of Eq. ([3| can be 
extended to include cavity damping and driving, as well 
as atomic damping via free space scattering of light (i.e., 
not into the cavity mode) as follows: 



d(c) 
dt 

d(a) 
dt 




(5) 



These equations include the additional non-unitary 
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damping and drive terms shown schematically in Fig. 
^b). The additional terms can be derived using input- 
output theory [52] and will now be discussed. But first, 
we note that at large fi-j-, the finite damping terms intro- 
duce small corrections to w± of order (T — k) 2 /J7^ [53] 
that may be safely neglected in this paper. 



1. Cavity damping and input- output fields 

As shown in Fig. [3](b) , the damping of the cavity field 
at rate k/2 — (ki + k 2 + «l)/2 is set by the mirror power 
transmission coefficients such that K\ % % — x /fsr- 
The total round-trip scattering and absorption fractional 
power losses at the mirrors L can be modelled by an 
additional beam splitter with field decay rate Kl = L x 
/fsr- Here, the cavity free spectral range is /fsr = c /2^, 
with 21 being the round-trip cavity length, and c the 
speed of light. 

The cavity field is driven by an external incident probe 
field described by the complex amplitude Cj. The re- 
flected and transmitted complex field amplitudes, c r and 
Ct respectively, will be detected to infer the number of 
atoms in The external field normalizations are cho- 
sen such that |ci irjt | 2 is the flux of incident, reflected, and 
transmitted probe photons in units of photons/second. 
The average number of incident, reflected, and transmit- 
ted photons Mi tTt t in a measurement time interval from 
t' to t' + T m is then 



Mi 



.M")fdt". 



(6) 



It is often convenient to express the number of probe 
photons coupled into the atoms-cavity system in terms 
of the measured "missing" photons in the reflected mode 
compared to the incident beam M m = Mi — M r . 

The reflected and transmitted fields can be found by 
first solving the coupled-driven Eq. ([5| for (c) and then 
using the results in the approximate relationships 

Cr = y/Kl (c) - Ci , 



Ct = (c) , (7) 

that hold in the limit of a high finesse cavity <C 1. 

2. Atomic damping via free space decay 

The atomic damping via scattering of light into free 
space (i.e., not into the cavity mode) is described by an 
effective amplitude damping rate T/2. To good approxi- 
mation, the probability decay rate T is simply the single- 
particle excited state |e) decay rate in free space [51] . The 
rate of scattering into free space is described by the field 
amplitude a s = y/T (a) , normalized such that the rate of 
photons scattered into free space is simply M s = \a s \ . 



The above picture of atomic damping can be further 
refined as shown in Fig. [3jb) . While the decay of exci- 
tation from the cavity mirrors is single-mode in nature, 
the atoms scatter light into many free space modes. This 
multimode scattering can be envisioned by replacing the 
single decay process via a single mirror with a weak beam 
splitter for each atom in [[). If the ensemble is optically 
thin along all directions except the cavity mode, then one 
can approximate that each atom decays into its own bath 



of states with an amplitude a s 



The 



normalization is necessary to impose the fact that if a sin- 
gle photon is emitted into free space, the probability that 
it was the ith atom that scattered the photon is reduced 
as the number of atoms increases. The total scatter- 
ing rate is the incoherent sum of the decay rates, repro- 
ducing the previous decay rate M s = T\(a)\ 2 . However, 
this refinement importantly emphasizes that the multi- 
mode free space scattering leads to in-principle informa- 
tion gain as to which particular atoms are in |f) , causing 
single particle collapse of the atomic wavefunction from 
a coherent superposition into an energy eigenstate, for 
example ( \\) + / v2 — > | f ) , thus destroying coher- 
ence. In contrast, the decay of light through the cavity 
mirrors leads to only collective information as to how 
many atoms total are in spin up and therefore preserves 
coherence. Thus information gained through the cavity 
will be useful for preparing conditionally spin-squeezed 
states, while the free space scattering is a competing de- 
coherence mechanism that serves to reduce the attainable 
degree of spin squeezing. 



3. Dressed mode linewidth 

The atomic and cavity damping cause the normal 
modes to decay, resulting in a mode broadening with 
approximate Lorentzian full width at half maximum 
(FWHM) resonance linewidth given by a weighted av- 
erage linewidth as 



n 2 t r 



4<4 



(8) 



For atom number counting, we find that the dressed cav- 
ity mode that is farthest from atomic resonance is most 
useful. For brevity, we refer to this mode's linewidth as 
simply k' such that k! — k' ± when \u>± \ > \u>^\- 



C. Full complex field response to probing 

The reflected and transmitted fields relative to the in- 
cident field c r j/ci = I rt t + iQ r ,t can be described in the 
complex plane by the real amplitudes I r>t and Q r ,t- We 
consider a single incident probe frequency close to reso- 
nance such that the probe frequency lo p is detuned from 
one of the two resonances by 8 p ± = (lo p —uj±). We assume 
the probe is near resonance 1 5 P ± \ <C uj + — u- such that 
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(a)Q 




(b) I 2 + Q 2 
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FIG. 5. (Color online) Transmitted and reflected probe elec- 
tric fields from driving the atoms-cavity system through the 
cavity mode, (a) The electric fields trace out circles in the 
I—Q phase plane as the probe detuning 5 P ± from the dressed 
cavity resonance varies from far-below to far-above resonance. 
The normalization is chosen such that the reflected electric 
field goes to 1 when far off resonance. The quantum noise 
of the probe normalized to the incident electric field is rep- 
resented as a fuzzy blob with rms diameter l/|c»|. In this 
illustration, a symmetric cavity fti = K2 is assumed, so that 
the circles have the same radius, (b) Corresponding power 
transmission and reflection. 



interference effects between normal modes can be ignored 
and that the modes are well resolved k' ± <C u>+ — w_ . The 
normalized transmitted electric field through the cavity 
is then to good approximation: 



l + (2S p± /n' ± Y 
_ f3 ± (2S p± / K > ± ) 
1 + (26 p± /k' ± ) 2 
The unitless amplitude (3± is given by 

2^/k 1 k 2 



(9) 



(10) 



(11) 



From Eq. |7]), the reflected field c r is just the sum of 
the transmitted field (rescaled for relative transmission 
coefficients) and the largely reflected field such that = 
1 - \] k\/k 2 1 u and Q r = ^Ki/^Qt- 

As shown in Fig. [5]ja), the phasor ct traces out a circle 
of radius f3±/2 in the complex plane as S p ± varies from -C 
k' ± /2 to ^> k'±/2. The translation I[ = I t — j3±/2 centers 
the circle traced out by the phasor I' t + %Q t at the origin. 
One then sees that the angle with respect to the real axis 
is given by ip± = arctan(Q t // t ') = arctan(2<5 p -|-/K' ± ). A 
similar translation holds for the reflected field. 



D. Probe vacuum noise and measurement 
resolution 

The size of the quantum vacuum noise that contributes 
uncertainty to measuring the position of the phasor is not 
changed by a linear transformation of coordinates in the 
complex plane. For our purposes, the noise can be de- 
scribed as a Gaussian probability distribution with equal 
and uncorrelated real and imaginary rms fluctuations of 
magnitude a v = 1/2. The rms quantum vacuum un- 
certainty Atp on the angle ip is then independent of the 
average value ip and is set only by the average number of 
detected photons in transmission Md — qdM t as 



AiP 



(12) 



The detection quantum efficiency qd includes any light 
loss and any excess technical or thermal noise of the 
detector relative to vacuum noise. The uncertainty 
Atp maps onto an uncertainty on the estimation of S p ± 
through AS p ± = \dS p ±/di/j\ Aip = K' ± Aip/2rj d . The de- 
tection sensitivity rjd is given by 



1 



1 + (2<S p± / K ' d 



(13) 



±) 



Probing near resonance S p ± = 0, one finds rjd = 1. For 
side-of- fringe probing S p ± — n'±/2, one finds rjd — 1/2. If 
the probe frequency is linearly and adiabatically scanned 
from S p ± <C k' ± to 8 p ± k' ± such that the total number 
of detected photons is fixed to the same Md as in the two 
previous scenarios, one finds rjd = 1/2. The optimal read- 
out assumes that as 6 P ± is changed, an adaptive homo- 
dyne readout is employed to maximize the measurement 
sensitivity to small changes in ip. In Ref. [TJ, heterodyne 
detection is employed so that adaptive detection is not re- 
quired. However, the effective quantum efficiency qd was 
reduced by 1/2 as a result of the heterodyne detection. 

It is straightforward to extend the analysis to a probe 
signal detected in reflection. However, one must parame- 
terize in terms of the measurable average number of miss- 
ing photons in the reflection port M m and the average 
number of incident photons Mi such that in Eq. ( 12 1, one 
substitutes M d -> {k 2 / Ki)M t q d {l + \jl- M m /Mi) 2 . 



III. QUANTUM-LIMITED SIGNAL-TO-NOISE 
AND FREE SPACE SCATTERING 

The measurement of the atomic population in |^) 
is achieved by precisely measuring the dressed mode fre- 
quency u+ or cj_ or some combination of the two. In 
essence, the approach used here converts the problem of 
measuring an atomic population into a frequency mea- 
surement. For atoms in a coherent superposition of |^) 
and quantum projection noise in the atomic popu- 
lation N+ causes the dressed mode frequency to fluctu- 
ate. In this section, we first derive the fluctuations on 
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FIG. 6. (Color online) Theoretical scaling of key quantities 
with cavity detuning 8 C expressed in units of the the collective 
vacuum Rabi frequency flf. (a) The dressed cavity linewidth 
k in units of the atomic excited state linewidth F. In (a) 
and (b), the scaling for different cavity finesses is shown for 
ft/r = 0.01, 1, and 100 (blue, red, and green, respectively), 
(b, right, black curve) The rms fluctuation of the dressed 
cavity mode frequency due to projection noise fluctuations 
Ao;^ roj decreases as 1/8 C above <5 c /f2f = 1. The normaliza- 
tion is chosen such that one should multiply by g/2^/2. (b, 
left, red, blue, green curves). The ratio of the projection noise 
fluctuation of the cavity mode to the dressed cavity linewidth 
Au} pio1 /k' is shown normalized such that the plotted values 
should be multiplied by g/(2\/2r). A large ratio is desir- 
able because technical noise may limit the ability to split the 
probed resonance by more than a fractional amount. 



the dressed mode frequency due to quantum projection 
noise as a function of cavity detuning S c . We then use 
the results of Sec. [TT] to obtain the average number of 
free-space-scattered photons per atom m™, when the 
measurement resolution on the probe field is sufficient to 
resolve the projection noise fluctuations of the mode fre- 
quency Aw pr °j. The quantity mP r °j is the key figure of 
merit that characterizes the degree to which a measure- 
ment is non-demolition. Three limits of cavity probing 
are identified, and a summary table of various key quan- 
tities in different regimes is presented. 



A. Projection-noise-driven fluctuations of mode 
frequencies 

As stated earlier, the atom number N-f can be de- 
termined by precisely measuring the dressed mode fre- 
quency uj + or u;_ with a cavity probe. The collective 



FIG. 7. (Color online) Theoretical scaling of key quantities 
with cavity detuning S c expressed in units of the the collective 
vacuum Rabi frequency Qf . In (a) and (b), the scaling for dif- 
ferent cavity finesses is shown giving k/F = 0.01, 1, and 100 
(blue, red, and green, respectively), (a, left, blue, red, green ) 
The average number of detected photons needed to resolve the 
projection noise fluctuations M pros normalized such that the 
plotted values should be multiplied by F 2 /(rjdg 2 )- (a, right, 
black) The ratio of the number of free-space-scattered pho- 
tons for every detected photon R s , normalized such that the 
plotted values should be multiplied by F /(qdn). (b) The cru- 
cial average number of scattered photons into free space per 
atom m^ ro1 when the atomic population measurement preci- 
sion is equal to the projection noise fluctuations. The normal- 
ization is such that the plotted values should be multiplied 
by r 2 /(16qN^g 2 ). In the bad-cavity limit of k » T (green 
curves), there is little fundamental advantage to operating 
away from resonance S c — 0. The technical requirements are 
simply increased as a result of detuning. As the finesse of the 
cavity F is increased, the amount of free space scattering falls 
roughly as 1/F until the good cavity regime is reached when 
k <C r (blue curves). Here, one must detune by roughly the 
critical detuning 8° in order to realize the full advantage of 
having increased the cavity finesse. Importantly, note that 
rris roj does not significantly decrease above 8° owing to can- 
cellation in this regime of the scaling of R s ~ 1 /8 2 with the 
scaling ofM d proi ~ 8 2 c . 



enhancement of the Rabi splitting by WNf produces a 
key enhancement of the measurement sensitivity needed 
to resolve projection noise. To concretely analyze the sig- 
nal to noise of the probing, we consider a measurement 
procedure most relevant to spectroscopy: we assume that 
for each experimental trial, all of the N ^> 1 total atoms 
arc initially prepared in spin down via optical pumping 
or otherwise. Each atom is then rotated into an equal su- 
perposition of spin up and down, preparing the ensemble 
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width at half maximum (HWHM) y^2g. Figure Mb) 
shows this scaling with detuning (black, left curve), xhe 
technical requirements on the experiment for resolving 
Au proj are increased with detuning. Other experimental 
imprecision and inaccuracies scale relative to the mode 
linewidth k' that one must split to the level of Aw pro -', 
therefore the ratio Aui pr ° s /k' is shown in Fig.[6|b, left) for 
three different bare cavity linewidths k/T — 0.01, 1, 100 
(blue, red, green). Note that in the good cavity limit 
s/r « 1 (blue), the experimental requirement on split- 
ting the mode line can be somewhat reduced at larger 
detuning owing to the rapid fall off of k' as 1/(5;? in the 
approximate region S c /Q^ G {1, 10}. 



FIG. 8. (Color online) The crucial average number of scat- 
tered photons into free space per atom such that the detected 
probe light allows resolution of projection noise fluctuations 
m^ r ° 3 versus the ratio of the cavity power to atomic pop- 
ulation decay rates k/T. The normalization is chosen such 
that the plotted values should be multiplied by T 2 / (16qN^g 2 ). 
Here, we assume that k is varied by changing the cavity finesse 
while holding the cavity length fixed (such that g is constant). 
Each trace represents a fixed ratio of the bare cavity detuning 
to the collective vacuum Rabi frequency 8 C /Qf , with values 
labeled on the traces and blue denoting S c /Qf > 0, green 
<5 c /fi-|- = and dashed yellow 5 C /Q.^ < . Here we consider 
probing the dressed mode at frequency u)+. At a fixed de- 
tuning (and therefore fixed projection-noise-driven frequency 
fluctuation size) a minimum is reached below which a better 
cavity is detrimental due to the dressed cavity linewidth k' 
becoming clamped while the ratio of free space scattering to 
detected photon number R a continues to rise. In the units 
above, the locus of the minimum is just 4k /T (red line). In- 
cluding the normalization, the locus of the minimum reduces 
to m5 roi = 1/ (qNfC). When including the dependence of g 2 , 
K ~ l/l on cavity length I, the minimum value of raf 01 is not 
changed by just shortening the cavity at fixed finesse. Only 
shortening the cavity length while simultaneously increasing 
the cavity finesse (such that k is constant) leads to a net fun- 
damental reduction of mJ roj . Finally, the vertical black line 
indicates the cavity linewidth for a finesse F = 10 6 cavity, 
near the highest currently achievable, for 87 Rb and I = 2 cm. 



in a coherent spin state (CSS). The population in spin up 
and down fluctuate about the average — = N/2 
with equal magnitude but perfectly anti-correlated pro- 
jection noise fluctuations ANf = — AiVj, = yN /2. 

The rms fluctuation Aw pr °j of the individual mode fre- 
quencies uj± caused by the projection-noise-driven fluc- 
tuations in iVf is found by linear expansion as Alu p1 '° 3 = 
(duj±/dN^) |jv t AN-\-. Evaluating using Eq. (fil), one finds 



Aw proj 



(14) 



The fluctuations of the two mode frequencies are equal 
in magnitude but opposite in sign such that the rms dif- 
ferential fluctuation is A(w + - cj_) proj = 2Aw proj . 

The projection noise variance (Aw pr °j) 2 decreases as a 
Lorentzian versus the bare cavity detuning S c with half 



B. Fundamental measurement noise and free space 
scattering at arbitrary detuning 5 C 

The resonance frequency uj of the farther detuned of 
the two dressed modes lj + or cj_ is measured relative to 
the known frequency of a coherent (and unsqueezed) laser 
probe. The rms uncertainty on the probe detuning Ai5 p 
is equal to the projection noise fluctuation level Alu p1 '° 3 
at an average detected photon number of 



Mr - 




(15) 



The passage of light through the cavity also leads to the 



scattering of M s 



T m probe photons into free space 



modes by the atoms in spin up. The ratio of free space 
scattered to detected photons R s = M B /M d is given by 
a weighted ratio of the two damping rates as 



1 



n 2 t 



r 

k 4w 2 



(16) 



The factor rj s plays an equivalent role to a quantum ef- 
ficiency and separately accounts for photons exiting the 
cavity via an undetected port. In the example we con- 
sider here, only the transmission port 2 is measured, and 
T] s = k 2 /k [see Fig. [3[ 



From Eqs. ( 16 ) and ( 15 ), one finds that the number of 
scattered photons into free space normalized to the total 
number of atoms N is 



R s Mr 3 



(17) 
(18) 



The single particle cooperativity parameter is C = 
(2g) 2 / kT and the total effective quantum efficiency is 




(19) 



For any arbitrary measurement resolution A5 P — 
aAw™ relative to the projection noise level, the re- 
quired average number of detected photons is simply 
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Aid = M^ 103 /a 2 , and the average number of scattered 
photons normalized to the total atom number N is m s = 
mf^/a 2 . 

A key result is that the free space scattering at a con- 
stant measurement precision a saturates to a finite value 
l/(4qN^Ca 2 ) as \6 C \ — > oo. This saturation 



is because in Eq. ( |18[ ), the ratio of free space to de- 
tected photons asymptotically decreases as 1/6%, but the 
required number of detected photons increases asymptot- 
ically as 6 2 . The non-demolition character of the mea- 
surement is ultimately set by the collective cooperativity 
parameter and quantum efficiency qN^C. This quantity 
physically sets the maximum rate at which collective in- 
formation can be extracted from the ensemble compared 
to the rate at which single-particle information is gained 
by the environment via multimode scattering of light into 
the many modes of free space. 

In the good cavity limit n <^ T, the frequency depen- 
dence 6 C of Eq. ( 18 ) can be understood in three regimes: 
the far-detuned dispersive regime |<5 C | > 6°, the near- 
detuned dispersive regime |<5 C | < 6°, and the resonant 
regime 6 C — 0. The critical cavity detuning 6° is given 
by 



6° 



rn t i 



IN+C 



(20) 



The critical detuning is the cavity detuning at which the 
dressed cavity lincwidth is k' = 2k, possible only in the 
good cavity limit k <C T. Expressions for the number 
of photons scattered into free space per atom mP roj , the 
absolute size of the projection noise fluctuations of the 
mode frequency Aw proj , and the dressed cavity linewidth 
k! are summarized in these different regimes in Table [I] 
Again, the quantity m^ ro) is critical for understanding 
the fundamental limits on both probe-induced heating of 
the sample and potential improvements on measurement 
sensitivity beyond the SQL. 



C. Minimizing at fixed maximum detuning 

In some experimental situations, a maximum probe de- 
tuning \6 C \ < (5 max is set by the energy structure of the 
atom. For instance, the ground state hyperfme splitting 
in 87 Rb imposes <5 max ps 6.8/2 GHz [7 . An optimum 
value of N^C can be found that minimizes m§ roj when 
\6 C \ = <5max- The scaling for this case is shown in the last 
line of Table [Tj Physically, the optimum value of NfC 
is reached (at a fixed detuning) when the dressed cav- 
ity linewidth is related to the bare cavity and atomic 
linewidths by k' /k — 2r/(r + n). In the resonant limit, 
6 C = 0, one finds an optimum k'/k = 1, while in the 
detuned limit one finds k 1 ~ 2k, i.e., the detuned cavity 
resonance is broadened by a factor of 2 at optimum. In 
this same limit, the ratio of rms fluctuation size to dressed 
cavity HWHM is given by Auj^/in' /2) = VC/2^2. 



IV. EXPERIMENTAL PROBING IN THE 
RESONANT LIMIT, 5 C = 

Here we consider the special case of probing in the res- 
onant limit, 6 C = 0, used in the experiment of Ref. [T]. 
In Sec. |IV A[ we give a brief review of the experimental 
setup. In Sec. IV B we consider the theoretical results 
of Sec. [III| in the limit of 6 C = 0. Sec. [IVC] and [TVD] de- 
tail technical advantages and sources of common-mode 
noise rejection obtained by working in the resonant limit 
and simultaneously probing both dressed cavity modes. 
Sec. |IV E| considers the effects of spatial mode mismatch 
of the probe modes in a standing-wave cavity on the 
correlation between the dressed mode frequencies. Fi- 



nally, in Sec. IV F we provide detailed descriptions of 
two schemes for probing the dressed mode frequencies at 
low probe powers. 



A. Review of experimental setup 

The two-level system whose populations we wish to 
measure are the clock states | f) = \F — 2,mp = 0) 
and \l) — \F — l,m F — 0) of the 87 Rb ground elec- 
tronic state 5 2 5 1 /2 with microwave transition frequency 
Whf/27r = 6.834 GHz [see Fig. |TJb)] . An ensemble of ap- 
proximately 10 6 87 Rb atoms are loaded into the TEMoo 
mode of an optical cavity whose frequency is made reso- 
nant with the |t) to |e) optical transition with wavelength 
A = 795 nm. Here, the electronic excited state is |e) = 
\5 2 Pi/ 2 ,F' = l,m F — 0). For 7r-polarized probe light, 
the states |f) and \e) form a nearly ideal two- level cavity- 
QED system because the matrix element between |f) and 
the only other nearby state \5 2 Pi/ 2 ,F — 2,mp — 0) is 
zero. Confinement and a fixed spatial registration of the 
atoms with respect to the cavity mode is accomplished 
via a 360 /iK deep intracavity 1-D optical lattice with 
wavelength Xi at — 823 nm. The radial temperature of the 
atoms is measured to be about 25 [iK. The rms extent of 
the atomic motion in the lattice is a: rms = y rms ~ 10 /im 
in the radial directions. In the axial direction, the atoms 
are trapped in the Lamb-Dicke regime with z rms 24 nm 
and Lamb-Dicke parameter w 0.2. 

The spatially averaged cavity cooperativity parame- 
ter is C = 3.2 x 10~ 3 , the cavity power decay rate 
is k — 2n(%.2 MHz), and the excited state |e) popu- 
lation decay rate is T = 27r(5.75 MHz). Over time, 
we found that the cavity linewidth has increased from 
9.2 MHz during the heterodyne Pound-Drever-Hall de- 
tection scheme described in Sec. IIVF II to 11.1 MHz 
during the heterodyne IQ detection scheme described 
in Sec. |IVF 2| The peak on-axis value of the coupling 

constant is 2g a = \J ^e}v m = 2tt(0.607 MHz). The 
atomic dipole moment is d ea — 1.465 C-m, c is the speed 
of light, e is the electric permittivity of vacuum, and 
h is Planck's constant. The effective mode volume for 
our symmetric standing- wave cavity is V m = -^ttwqI = 
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TABLE I. Regimes of cavity probing. The regime name and assumptions used to define the regime are provided in the first 
and last columns. The quantity m^ 10 ' is the average number of photons scattered into free space normalized to the total atom 
number N, required to resolve an rms fluctuation *jNf/2 in the spin up population equal to the projection noise level. The 
quantity Aw proj is the rms angular frequency fluctuation of a single coupled atoms-cavity mode ui± due to projection noise. 
The quantity «' is the dressed cavity power decay linewidth, here taken for the mode detuned farthest from atomic resonance. 
The single-particle cooperativity C, the number of atoms in spin up Nf = N/2, the single-particle cavity coupling g, the empty 
cavity power and atomic population decay rates re and V respectively, and the collective vacuum Rabi frequency fit are related 
by the following: fit = U Nf2g, NfC = Nf(2g) 2 / ' nT . The detuning of the empty cavity resonance frequency from the atomic 
transition frequency is 8 C , and the critical detuning at which re' = 2re is (<5°) 2 = fi 2 F/4re, assuming the good cavity limit F ^> re. 
The maximally detuned regime assumes that the quantity NfC is chosen to minimize m^ T01 in the presence of the constraint 
that the cavity detuning cannot be made larger than some maximum value <5 max set by technical constraints on resolving the 
projection noise fluctuations or fundamental constraints set by the internal energy level structure of the atoms being probed 
(for instance the ground state hyperfine splitting in 87 Rb). 
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0.0748 mm 3 . The mirror separation is I = 1.91 cm, the 
mode waist wq = 70.6 /im, and the mode Rayleigh range 
is zr = 1.97 cm. The axially averaged coupling constant 
is 2g = 27r(0.526 MHz). Because of the tight confinement 
of the atoms compared to the waist of the probe mode 
and the periodicity of the standing-wave probe mode, 
thermal averaging only slightly reduces 2g by 4%. 



B. Fundamental scalings related to projection noise 

On resonance 5 C = 0, the absolute size of the projection 
noise fluctuations is maximized, i.e., 

AtJ p™; = 9 (21) 

Note that the rms fluctuation is independent of N, as 
experimentally demonstrated in Fig. [9j The same is true 
for the FWHM linewidth which is simply equal to the av- 
erage linewidth [53] due to the equal photonic and atomic 
contributions to the normal modes: 

«4 = (K + r)/2. (22) 

To be able to resolve projection noise, one must detect, 
on average, a number of probe photons in transmission 



given by 



The ratio of free-space-scattered photons to detected 
probe photons in transmission is 

1 r 

R s ± = ■ (24) 

QdVs K 

Finally, the number of scattered photons into free space 
for measurement uncertainty at the projection noise level 
normalized to the total number of atoms N is 



1 / r x 2 

proj 1 I - L 



2qNC 



(25) 



If the cavity length and mode volume are fixed by ex- 
perimental constraints, then one is, in principle, free to 



minimize Eq. ( 25 ) by varying the finesse of the cavity 
mirrors, until a minimum value of mP roj = 2/qNC is 
reached when k = T. The minimization with respect to 
cavity finesse accounts for the fact that the cooperativity 
C scales as 1/re. 
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FIG. 9. (Color online) In the resonant limit S c = 0, the 
rms projection noise fluctuation of the population Nf is inde- 
pendent of total atom number N. Shown here are the mea- 



(i) 



sured rms fluctuations of the dressed cavity modes A(fii 

Q.^)/2ty after subtracting out the probe measurement noise 
in quadrature. The data agrees with the a priori prediction 
(line) for the projection noise level s/2g/2n = 0.358 MHz 
very well. At larger atom numbers N > 7 x 10 5 , additional 
technical noise was observed. In contrast, probing in the off- 
resonant regime results in a cavity mode frequency noise vari- 
ance that scales linearly with N, and is always smaller than in 
the resonant limit [see Eq. |14l]. To explore potential added 
noise from the 7r-pulse rotation that swaps the two atomic 
populations Nf and N± between the two splitting measure- 
ments, we apply the rotation about an axis perpendicular to 
the Bloch vector (black points), and about an axis parallel to 
the Bloch vector (red diamonds). 



C. Noise rejection in the resonant limit 



Section [IVB] has outlined fundamental quantum limits 
on detection in the S c = regime. In this subsection, we 
describe how the S c = regime is useful for rejecting 
other sources of technical noise. Simultaneous probing 
of both transmission resonances ui± is advantageous for 
canceling technical, but important, laser frequency noise 
and cavity frequency noise. Importantly, the previous 
analysis for the value of m^ mi is not modified if M^ roj /2 
photons are used to measure each mode frequency uj + 
and uj- . 

From Eq. Q, we find that the measured differential 



quantity uj+ 



a 



sensitive to the cavity detuning S c 



is only quadratically 
Note that the abso- 



lute sensitivity to 5 C is decreasing as l/\ 



a favorable 



scaling for working with large atom number. As an ex- 
ample, a typical Q^/2ir = 300 MHz and a detuning of 
5 c /2ir = 1 MHz causes a change in the measured split- 
ting 100 times smaller than projection noise fluctuation 
in the collective Rabi frequency Af^ roj . 

Simultaneous or near-simultaneous measurements of 
both resonances can strongly suppress the impact of laser 



frequency noise. To make these measurements, one can 
use low phase noise RF modulation of a single laser with 
much higher phase noise to generate probe components 
near w+ and w_. The laser carrier frequency noise is 
common-mode to the two probes and cancels in the dif- 
ference uj + — lu- . The measured power spectral density 
(PSD) of the probe-cavity relative frequency noise S v (f), 
shown in Fig. |10[ was measured for our system using 
side-of-fringe transmission. If both the upper and lower 
transmission resonances are measured for time T m with 
time separation r m between the measurements, the noise 
variance in the measured splitting due to laser noise is 



/ AD! r asor ' 



2?r 



S v {f)T{f)df. 



where the transfer function T(f) is 

4sin 2 (7r/r m )sin 2 (7r/T m ) 



T(f) 



(26) 



(27) 



= Av/n. Integrating Eq 
= S^Tm/T^ for good time 



(|26|), we 

over- 



If the laser possesses a Lorentzian lineshape with 
FWHM Av in Hz as measured in heterodyne ver- 
sus a perfect reference laser, the PSD is constant 
with value S u (f) 
find (Af^ asor /27r) 2 
lap r m < T m , and S v /T m for no time overlap r m > 
T m . Completely unoverlapped measurements would yield 
an absolute measurement uncertainty of Af2| ascr /27r = 

U Ai>B m /iT, the geometric mean of the measurement 
bandwidth B m = l/T m and the laser FWHM linewidth 
Az/. For a typical extended-cavity-diode-laser Av = 
100 kHz, and Afi| asor /2tt w 180 kH Zl /^s/T m . For partial 

time overlap Arj' ascr /27r 180 kH Z% /^s^T m /T^. 

In practice, lasers exhibit low frequency 1// noise, but 
so long as l/r m is much larger than the noise corner 
frequency, the above description is adequate. The laser 
noise can be reduced by spreading the Md detected probe 
photons over a longer measurement period, but at the 
possible expense of additional decoherence of the atomic 
quantum state and the potential additional contributions 
from 1// laser and cavity frequency noise. 



D. ac Stark shift cancellation 

The intracavity probe light induces an ac Stark or light 
shift of the |^) state. At resonance, simultaneous probing 
of both modes uj± leads to a cancellation of the shift as 
the induced shifts have opposite signs as the two modes 
are symmetrically spaced above and below the atomic 
transition frequency uj^ e . In the case of inhomogeneous 
coupling of the atoms to the cavity mode (as was the case 
in Ref. [1]), the ac Stark shifts are inhomogeneous and 
lead to a loss of signal and a reduction in the effective- 
ness of composite pulse sequences used for manipulating 
the ground state spin manifold. The integrated phase 
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FIG. 10. (Color online) (a) Measured probe laser frequency 
noise spectrum S v (f) (solid black curve) and for reference 
the calculated white frequency noise level for a Lorentzian 
50 kHz FWHM linewidth laser (horizontal black line). Also 
shown are theoretical transfer functions relating laser fre- 
quency noise to Rabi splitting noise for T m = 20/j.s and 
Tm/T m — 0.05,0.1, 1000 (dotted green, blue, and red curves, 
respectively) The inset of (b) defines the measurement time- 
ing. For the case T m /T m = 1000, the sin 2 (7r/T m )/(7r/T m ) 2 
envelope of the transfer function is shown here, obtained by 
averaging over the underlying fast oscillations, (b) RMS noise 
in the Rabi splitting noise Afi!|. aser /27r due to the measured 
laser frequency noise versus r m for T m = 20 fjs (solid blue 
curve) and for T m = 200 /us (solid red curve). For a white 
laser frequency noise spectrum of the level shown in (a), the 
integrated Rabi splitting noise spectrum versus r m is shown 
in dashed blue and red curves. Vertical dashed lines indicate 
the two different T m used for the transfer function integration. 
Most of the integrated noise at larger r m can be attributed 
to the laser frequency noise bump near 17 kHz. 



shift caused by passing probe light through the cavity 
was measured using a microwave Ramsey sequence and 
observing the phase shift of the Ramsey fringe versus the 
number of photons coupled into the atoms-cavity system. 
In Fig. [TTJ we show the probe-induced phase shift if only 
the w+ mode or the o>_ mode is probed. We then show 
the high degree of cancellation of the phase shift, and 
therefore the average induced ac Stark shift, when both 
the ui+ and cj_ modes are probed simultaneously. 



E. Spatial mode matching of atoms-cavity modes 

In a standing-wave cavity, the u± modes arise due 
to coupling of the cavity mode to atoms located away 
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FIG. 11. (Color online) The relative quantum phase between 
states |t) and \\) is light shifted by passing probe photons 
through the atoms-cavity system near only the upper normal 
mode frequency ui+ (blue solid circles) or only near the lower 
normal mode frequency u_ (red open circles) . The shifts have 
opposite sign and cancel if probe photons are simultaneously 
passed through the atoms-cavity system at both frequencies 
(purple squares). The measurements are performed using a 
variable number of photons coupled into the atoms-cavity sys- 
tem M m during a microwave spin-echo sequence (inset) em- 
ployed to remove background dephasing from the trapping 
optical lattice. 



from the nodes of the standing-wave of each mode. By 
probing at two different frequencies, it is possible that 
the two standing-wave probe modes effectively couple to 
two slightly different sub-ensembles of atoms. Projection 
noise fluctuations in the two sub-ensembles are not per- 
fectly correlated such that the fluctuations of w+ and w_ 
are also no longer perfectly correlated. The contribution 
of projection noise to the difference of the mode frequen- 
cies w+ — ui- would then be reduced. In the extreme 
limit of the nodes of the lu + mode being aligned with the 
antinodes of the w_ mode, the two sub-ensembles are 
completely independent, the projection noise driven rms 
fluctuations of the difference A(uj + — cj_) would be v2 
lower than for the case of perfect overlap. Here we calcu- 
late this effect and provide numerical results indicating 
that this effect is small in the experiment of Ref . [T] . 

In vacuum, the wave vectors for resonant probe light 
k± = lj±/c have a differential wave vector given by 
k + — = Q^/c. However, if the atoms uniformly fill 
the cavity mode along the axial or z direction, then the 
index of refraction of the atoms at the resonant frequen- 
cies oj± must be such that two wavevectors k' ± = n±k± 
are identical inside the atomic medium. This is guaran- 
teed by the round trip phase being a multiple integer of 
27T at resonance. 

If, however, the atoms do not uniformly fill the cav- 
ity, then the wavevectors of the probe modes in vacuum 
and inside the atomic ensemble can be different, leading 
to coupling to slightly different atomic sub-ensembles. 
Using the original condition, one finds the indices of re- 
fraction of the atomic ensemble at position z along the 
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FIG. 12. (Color online) Illustration of the concept of spatial 
mode matching of the atoms-cavity modes, (a) In vacuum, 
the electric field of the higher frequency mode u+ (blue) has a 
shorter wavelength than the lower frequency mode oj_ (red), 
leading to a phase advance of the u+ mode relative to the 
oj- mode as it advances from the left mirror to the left edge 
of the atomic medium (green). In the atomic medium, the 
refractive index n± at frequency u± is such that the co+ mode 
has a longer wavelength than the ui- mode so that the two 
modes phase up in the middle of the cavity. For reference, 
the electric field for the bare cavity mode uj c (dashed black) is 
shown. The phase mismatch in this illustration is deliberately 
made large by choosing — /fsr- The length of the 

uniformly distributed atomic medium is set to d = 1/4, again 
for illustration, (b) The spatial phase deviation <f>± of the ui± 
modes from the empty cavity spatial phase. The red(blue) 
curves corresponding to <j>+ and <j>- respectively are used to 
generate the standing- wave modes in (a). 



In the experiment of Ref. [T] , the atoms were roughly 
Gaussian distributed along the cavity axis with standard 
deviation a p = 0.84(9) mm, much less than the cavity 
length I = 1.91 cm. In this case, the spatial phase devi- 
ation of the standing-wave probe modes from the empty 
cavity spatial phase in the atomic medium is given by 



4>±{z) = ± 




(30) 



where Erf(-) is the error function and \z\ < d/2. 

At z = CT p , and typical Q t = 2tt(300 MHz), the differ- 
ential probe spatial phase 4>+{z) — 4>-{z) is expected to 
be only 31(3) mrad. Numerical simulations of the phase 
mismatch and projection noise show that this phase mis- 
match leads to < 1% reduction in the variance of the 
difference lu + — o;_ relative to the variance in the same 
quantity for perfect spatial overlap of the two modes with 
one another. 



F. Simultaneous mode probing at low light powers 

We implemented two schemes to simultaneously mea- 
sure the two modes of the vacuum Rabi splitting. 
The first is based on a modified Pound-Drever-Hall 
(PDH) scheme employing heterodyne detection to pro- 
vide quantum- limited readout. The second technique is 
a heterodyne probe detection with full IQ quadrature 
measurement combined with active servoing of the opti- 
cal path length. 



cavity axis evaluated for probe frequencies at the nominal 
resonance frequencies ui± are 



n± (z) = 1 T 



1 



pi z ) 



(28) 



4/fsr k± 

where /fsr = c/2l is the cavity free spectral range, and 
p-\ (z) is the coarse-grained linear probability distribution 
function that an atom in \f) is located at z and normal- 
ized such that integrating over all of space yields 1. The 
coarse graining is performed over many lattice sites, and 
quantum fluctuations are ignored in this expression as 
they are a higher order effect. 

As a first application, we assume the atoms uniformly 
fill a region of space of length d centered in the middle of 
the cavity of length I. In the atomic medium, the spatial 
phase deviation <f>± (z) of the standing- wave electric fields 
cos {oj c z I c + cf>± ) from the empty cavity spatial phase is 
given by 



zVL 



t 



1 



for 



d 

< _. 
~ 2 



(29) 



"4/fsr \l 

In the limit of uniform filling d = I, the two probe modes 
are spatially matched, otherwise they acquire a differen- 
tial spatial phase 4>+{z) — 4>—{z) that grows linearly as 
z. 



1. Heterodyne Pound-Drever-Hall detection 

The vacuum Rabi splitting is measured by sweeping 
the probe laser frequency and detecting the light reflected 
from the cavity input mirror. The resonance signals are 
extracted using a modified PDH detection scheme that 
employs optical heterodyne detection to achieve near 
photon-shot-noise-limited detection at probe powers of 



< 100 pW to several nW (see Fig. 13 ). To reduce sensitiv- 
ity to laser and cavity noise, we measure both resonances 
uj + and cj_ simultaneously as discussed in Sec.|IVC] The 



two probe components are generated by phase modulat- 
ing the probe light to create sidebands at ±151 MHz from 
the carrier. The modulation frequency is chosen to match 
the vacuum Rabi splitting of 302 MHz. Additional phase 
reference sidebands at ±15.3 MHz are then imposed by 
a second phase modulator. To avoid degeneracy of the 
two probe components in the RF spectrum generated at 
the photodiode, the probe carrier is blue-shifted 52 MHz 
with respect to the reference heterodyne beam derived 
from the same laser as the probe. 

The light reflected from the cavity is overlapped with 
the heterodyne light on a 90/10 beam splitter. Both 
beams are coupled into a single mode polarization- 
maintaining (PM) fiber to ensure high spatial mode 
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FIG. 13. (Color online) (a) Experimental setup of the het- 
erodyne Pound-Drever-Hall (PDH) scheme used to probe the 
vacuum Rabi splitting. Both components of the vacuum 
Rabi splitting are simultaneously probed using the first order 
phase modulation sidebands generated by the first electro- 
optic modulator (EOM) and spaced ±151 MHz from the car- 
rier. Additional ±15.3 MHz PDH phase reference sidebands 
are imposed by the second EOM. (b) Relevant energy level 
structure of 87 Rb. The bare cavity is made resonant with the 
t) to |e) transition, and 7r-polarized probe light is used, (c) 
Simultaneous probing of the vacuum Rabi splitting. The blue 
and red points are the individual PDH error signals simul- 
taneously measured by sweeping the phase-modulated probe 
laser frequency across the higher and lower frequency modes 
of the vacuum Rabi splitting respectively. The blue and red 
solid lines are curve fits to the data using an analytic model 
of the PDH error signal. 



matching. The heterodyne laser power delivered to the 
photodiode is typically 0.8 mW. Working at low het- 
erodyne powers provides additional immunity to clas- 
sical intensity noise, as well as increased probe detec- 
tion efficiency. The photodetector consists of a Hama- 
matsu S5973 photodiode, ac coupled to an Analog De- 
vices AD8015 250 MHz transimpedance amplifier. 

The two probe components are separated in the het- 
erodyne RF spectrum using high and low pass filters, 
and notch filters at 99 MHz and 203 MHz. Each time- 
domain signal is squared using an Analog Devices AD835 
250 MHz voltage multiplier. This step essentially mimics 
the square law detection of the photodiode that occurs in 
the usual PDH detection scheme. The resulting 15.3 MHz 
signals are then demodulated to dc. Figure |l3|c) shows 
the resulting nearly simultaneously measured PDH error 
signals when the probe laser frequency is swept across a 
vacuum Rabi splitting. Linear fits to the central 15 /is of 
each probe component sweeping through resonance were 
used to extract the vacuum Rabi splitting fi-f . 

Using this technique at low probe powers, we were able 
to measure the collective spin at a precision of 3.5(4) dB 
above the quantum projection noise level for an ensemble 



of 7.1(3) x 10 5 atoms with a calculated m s = 0.16 pho- 
tons/atom worth of free space scattering, directly pro- 
portional to photon recoil heating and loss of contrast as 
described in Sec. |IVB| Simple extrapolation implies that 
with ss 2.2 x more probe photons, the measurement pre- 
cision would reach the projection noise level while causing 
~ 0.35 photon/atom worth of free space scattering. This 
PDH probing technique suffers from added noised dur- 
ing the electronic squaring process, degrading the signal- 
to-noise. However, rather than further optimizing this 
approach, we chose to pursue a more sensitive readout 
scheme described in Sec. IIV F~2l below. 



2. Heterodyne IQ detection 

For the results of Ref. pQ, we employed heterodyne 
detection as described in Sec. |IV F 1[ but with the probe 
sidebands directly IQ demodulated so that the full ampli- 
tude and phase response during a frequency sweep could 
be used to extract the mode resonance frequencies, as de- 
scribed in Sec. |II C| The probe beam was only phase mod- 
ulated at 150 MHz; otherwise the rf configuration was the 



same as in Sec. IV F 1 The phase modulation generates a 



high frequency and a low frequency probe sideband, blue 
and red colored beams respectively in Fig. 14]Ja). To 
help remove optical path length noise in the heterodyne 
interferometer, we also detected the probe carrier and 
stabilized its phase measured at the photodiode using a 
servo loop with the 106 MHz acousto-optic modulator as 
the phase shifter. 

We detect both the reflected and transmitted probe 
light and hence achieve better signal-to-noise in our 
symmetric cavity system as compared to just detecting 
only one port. Using IQ demodulators (Analog Devices 
ADL5387 for the high and low frequency probe side- 
bands, and Linear Technology LT5517 for the carrier), 
we phase coherently demodulate to obtain full ampli- 
tude and phase information of the probe electric field, 
or equivalently the full in-phase (I) and quadrature (Q) 
amplitudes. Figure 14 ^b) shows typical signals recorded 



in both reflection and transmission. The high (low) fre- 
quency probe component sweeps across the resonance ui + 
(cj_), with the resulting signals shown in blue (red) in 
Fig.^b, c, d). 

Phase stabilization is important for getting a good 
phase estimate for the IQ fit in reflection because the 
sweep spends relatively more time far from resonance 
than near resonance. Thus any phase wander of the het- 
erodyne interferometer will show up as a phase wander of 
the reflected electric field. On the other hand, the phase 
stabilization is not crucial for the IQ transmission fit be- 
cause there is no transmission off resonance, and there- 
fore the signal is in the detection noise and has no well- 
defined phase anyway. However, the time to sweep across 
resonance needs to be fast compared to the timescale for 
phase wander so that the circle defined by the IQ trans- 
mission response does not rotate due to phase drift as 
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efficiencies in reflection and transmission. With an av- 
erag e num ber of missing probe photons in reflection (see 
Sec. II D) M m — 0.96(5) x 10 5 probe photons for a sin- 
gle measurement of , we achieve a probe measurement 
precision 6.5(5) dB below the projection noise level for 
Nq = 7.0(3) x 10 5 atoms obtained from the relative noise 
between two repeated measurements on the same 
atomic ensemble. 

In connection to the theory presentation in Sec. |VI| 
we outline here the spectroscopic enhancement gained 
from the QND measurement in our experiment. The 
population difference and hence spin projection J z = 
(iV-j- — N^)/2 is determined from a measurement of Q-j-, 
followed by a coherent swapping of populations between 
spin up and down using a microwave 7r-pulse, after which 
f2^ is measured again. From repeated measurements 
of J z , we observed a conditional spin noise reduction 
in the difference of two successive QND measurements 



A 



{^J^ — jJ^J of 2.6(3) dB below the projection noise 
level. By subtracting out probe measurement noise vari- 
ance from the conditional spin noise variance, we in- 
ferred that a state with 4.9(6) dB reduced spin noise has 
been prepared. After the first measurement the 
measured contrast of the Ramsey fringe decreased from 
Ci = 0.97(1) before the measurement to C f = 0.82(2). 
Taking into account the 1.5 dB loss of contrast, we con- 
cluded that a conditionally spin-squeezed state with spec- 
troscopic enhancement 3.4(6) dB relative to the SQL has 
been prepared in an ensemble of 7.0(3) x 10 5 atoms . 



FIG. 14. (Color online) (a) Experimental setup of the het- 
erodyne IQ detection scheme used to probe the vacuum Rabi 
splitting. Both components of the vacuum Rabi splitting are 
simultaneously probed using the first order phase modulation 
sidebands spaced ±150 MHz from the carrier. Probe trans- 
mission is also detected using a similar setup as the reflec- 
tion detection (not shown), (b, c, d) Simultaneous probing 
of the vacuum Rabi splitting in transmission and reflection. 
The blue and red data points are the individual Rabi split- 
ting signals simultaneously measured by sweeping the phase- 
modulated probe laser frequency. The resonance frequencies 
are extracted from fits (black curves) to the data plotted para- 
metrically versus the probe frequency in the IQ plane (b). A 
rotation has been applied to align the circles along the I- 
quadrature for display purposes. The corresponding power 
I 2 + Q 2 (points - data, curve - from IQ fit results) (c) , and 
phase shift arctan(Q/I) (points - data, curve - from IQ fit 
results) for both transmission and reflection (d). 



V. QUANTUM BACK-ACTION LIMITS ON 
DETERMINING J z 

Having presented experimental examples of probing an 
atomic ensemble in the resonant limit near and below 
the quantum projection noise level in Sec. |IVF| we now 
study the limitations on measurement resolution arising 
from Raman spin flips caused by free space scattering. 
We will begin by considering a simple three level model 
that will be used in Sec. |VI| to calculate spectroscopic 
enhancements relative to the SQL. The simple model will 



be extended in Sec. VII to probing the clock and cycling 
transitions in 87 Rb. 



Average impact on J z due to free space 
scattering 



the probe sweeps across resonance. 

The resonance frequencies were extracted using least 
squares fitting in the IQ plane, yielding estimates fl T ^ efi 
(f2* rans ) of the Rabi splitting from the IQ reflection 
(transmission) fit, respectively. Empirically, we find that 
the weighted average (0.61f^ ofi + 0.39fij rans ) exhibited 
the least noise in two repeated measurements of the same 
Rabi splitting. This ratio reflects the relative quantum 



1. Simple three-level model 

In this subsection, we consider how the free space scat- 
tering changes atomic population in the spin up and 
spin down two-level manifold through Raman or spin 
flip events. We define collective spin operators J x ,y,z = 
52i=i 4y zi wn ere Jx,y,z ^ s the single-atom spin operator 
for the ith atom such that J* | ti) = § 1 ti) j J % z \ U) = — \ I U 
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), etc. The collective spin operator is J z = (Nf — N±)/2 
where and N± are the atomic population operators 
defined in Sec. |II A| Expectation values of the collective 

spin operators are denoted by J x ,y.z = (jx,y,z\ The 
collective Bloch vector is J = (J x , J y , J z ). The radius of 



the collective Bloch sphere is 1Z = 




. For arriving 



at the results presented in this section, only the atomic 
populations matter, and coherences are irrelevant. 

A detailed accounting of all possible Raman scattering 
processes is required to make accurate predictions for a 
multi-level atom. We consider here the simplest model 
that captures the essential physics. In this toy model, the 
only states in the problem are the two-level system |f) , \\) 
and the optical excited state |e), as described in Fig.[IJb). 
We assume that a free space scattering event causes an 
atom to spin flip from |f) to |4-) via the intermediate state 
|e) with probability p. It is easy to see that the average 
change in the spin up population per photon scattered 
into free space is dN^/dM s — —p, and correspondingly 
d.J z /dM s = -p. 



2. Experimental verification of the average impact on J z 
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FIG. 15. (Color online) Measurement of change in Rabi 
splitting (fil — Oi )/27r due to free space scattering ver- 
sus the missing photons M m in the reflected port. The 
slope is predicted to be ' Using the indepen- 

dently and accurately known values of g — 2n(253 kHz) and 
N t = 3.50(15) x 10 5 , we inferred M s /M m = 0.408(4), in 
excellent agreement with the theory prediction of 0.405(6). 
The inset shows the experimental sequence used to measure 
- flf\ A CSS state is prepared using the n/2 pulse 

(dashed line), and then the first measurement fij ' is made. 
Next, a variable number of probe pulses are applied to cause 
free space scattering, and M m is measured. Finally, a second 
measurement fii is made to determine Qi. — 



To connect with experiments performed in the reso- 
nant limit 5 C = 0, where the Rabi splitting = y / iV^2g 
is measured directly rather than J z , the previous result 
can be re-expressed as the change in the Rabi splitting 
per photon scattered as 



dft t 
dAL 



= ~P- 



(31) 



Accounting for the multi-level structure of 87 Rb and 
inhomogeneous coupling effects, we find in our system 



dCl-f 
dM s 



7 



9 



(32) 



It is more convenient to measure the "missing" probe 
photons from the reflected port M m than the multimode 
free-space-scattered photons M s . The quantities are re- 



lated by the following, using Eq. (Ill, 



B. Measurement Resolution on J z in Presence of 
Quantum Back-action 

1. Diffusion of J z 

Free space scattering causes J z to change on average 
by a certain amount, while the random nature of the 
spin flip process leads to a random walk or diffusion of 
the collective Bloch vector's spin projection J z . Provided 
multiple scattering can be neglected, i.e. pm s <C 1, the 
diffusion process can be described by the relation 



(J z (m s ) - J z (0)Y 



(AJ : 



z,CSS 



Apm s , 



(34) 



with the spin flip probability setting the diffusion con- 
stant Ap, and the random variable J z (m s ) describing 
the z-component of the Bloch vector after m s scatter- 
ing events per atom. The diffusion is normalized to the 
projection noise level for a CSS (AJ^css) 2 = N/4. 



r + k — Ki 



-Mr, 



(33) 



2. Measurement imprecision due to photon shot noise 



which evaluates to M s = 0.405(6)M m for our system pQ. 

Experimentally, we observed a linear decrease in the 
Rabi splitting versus the "missing" photons M m (see 
Fig. 15 1. From the slope d£l^-/dM m , we inferred M s — 
0.408[4jM m , confirming the rate equation analysis used 
in the derivation of Eq. (|32| . 



The measurement imprecision AJ™ cas is a weighted 
average of the state preparation variance or prior knowl- 
edge, and the probe vacuum noise: 



AJ" 



A, I 



z.CSS 
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(35) 
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FIG. 16. (Color online) Normalized spin noise variance ver- 
sus fractional free space scattering m s /m^ 101 . The curves as- 
sume no prior knowledge, i.e., Cprcp — > oo, and NC = 10 3 . 
The measurement variance decreases as l/m s from averaging 
down photon shot noise (dotted black line) until noise due 
to random spin flips (dashed red and black lines) takes over 



at large m a . The minimum noise variance is 4^/pm? roj 



\/9>p/NC assuming perfect quantum efficiency and probing 
in the far-detuned limit. This minimum is r eached wh en 
two noise contributions are equal at m° pt = J mf oi /4p — ► 

1 / ^/8pNC assuming perfect quantum efficiency and probing 
in the far-detuned limit. A larger collective cooperativity NC 
or lower spin flip probability p allows photon shot noise to be 
averaged down more before spin flip noise takes over, as il- 
lustrated by the two curves for p/NC = 1CP 5 (black) and 
1CP 3 (red). The locus of the minimum variance is 2ro a /mJ roi 
(green). Constructing an optimal Kalman filter to exponen- 
tially weight the measurement record at larger m 3 more than 
at smaller m s allows a factor of two improvement (K = 2) 
in the minimum noise variance to ^j2p/NC, a limit that is 
asymptotically approached. 



The first term arises from the ability to split the reso- 
nance line due to photon shot noise, which scales as mj 1 . 
The parameter Cprcp = (AJ ZjPrcp /A J Zi css) 2 character- 
izes the state preparation noise or prior knowledge rela- 
tive to the projection noise level. For example, Cprcp = 1 
corresponds to preparing a CSS with no added techni- 
cal noise, Cprcp - ► oo corresponds to preparing a state 
with essentially unknown J z prior to measurement, and 
Cprcp < 1 corresponds to a prepared state with reduced 
spin noise. 



J z or equivalently the atomic population Nf? As a first 
pass, the total variance (AJ Z ) 2 in the estimate of J z is 



found by adding the measurement imprecision [Eq. ( 35 )] 
and spin- flip-induced diffusion of J z [Eq. (34)]: 



AJ' 



AJ. 



z.CSS 



TT^s J Cprcp 



4pm s . (36) 



Averaging down photon shot noise determines the z- 
projection of the Bloch vector more and more precisely. 
Eventually, however, scattering-induced diffusion of J z 
causes the value of J z measured at earlier times to be- 
come less correlated with the value of J z measured at 
later ti mes , adding noise to the estimate of J z , as shown 
in Fig. 16 The optimal resolution AJ° pt occurs at an 



optimal scattering m° pt where the noise contributions 
from measurement imprecision and diffusion due to Ra- 
man spin flips are equal: 



m° pt = 



1 



1 



1 
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AJ 



y/8pqNCh(S c ) Cprc P 2qNCh(S c ) 
1 / 8p 1 
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(37) 
(38) 



where the constant K = 1 in the present analysis, and 
K = 2 in an optimal Kalman filter analysis considered in 
Sec. VB4 Unless otherwise stated, K = 1 is assumed 
for the rest of this paper including figures. The fac- 



tor ^4(Cprep,<$c 



€ [1, 2] is a small 



improvement due to including prior knowledge of the 
state. For Cprcp < \Z^p/qNCh(6 c ), the optimal scat- 
tering m° pt goes to zero because we already have such 
good prior knowledge that measuring the state does more 
harm than good through introducing Raman scattering- 
induced spin-flips. The detuning dependence has been 
lumped into 



h(S c ) = 



to(6 c ) 



(39) 



where k'(8 c ) is the dressed cavity mode linewidth in- 
troduced in Eq. Q, oj(S c ) is the dressed cavity mode 
frequency introduced in Eq. |4), and fl-f = y/N^2g is 
the collective vacuum Rabi splitting. The parameter 
h(S°) — 1/4 at the critical detuning, and h(S c ) — > 1 as 
\5 C \ — > oo. In the limit of no prior knowledge Cprcp — > oo, 
no Kalman optimal filtering K = 1, and probing in the 
far-detuned limit, Raman spin flip noise limits the achiev- 



able resolution to 



3. Balancing measurement imprecision against diffusion of 

Jz 

Given the diffusion in J z , and the photon shot noise 
in the probe, we must ask: How well can one measure 



4- Optimal Kalman estimator of J z 

More rigorously, one can construct a Kalman filter to 
provide a minimum-mean-square-error estimator for J z 
at the end of a continuous measurement record |231 1551 
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I5B] . The measurement record is best parametrized by 
the average of free-space-scattered photons per atom m s 
instead of time. Balancing between averaging down pho- 
ton shot noise and increasing decorrelation due to spin 
flips, the optimal Kalman filter is a simple exponential 
de-weighting of the measurement record at smaller m s 
relative to larger m s . The exponential weighting time 
constant tw is given by the amount of time tw it takes 

for an average rnf — yj m^° 3 /Ap — 1/ \/8pqNCh(6 c ) 
scattered photons per atom to occur; T\y can be obtained 
from Eq. (16). The assumptions that the length of the 



measurement record is much longer than tw and that 
<C 1 have been made. In the Fourier domain, the 
filter is equivalent to a single-pole low-pass filter with 
corner frequency % = 1 /tw . Implementing the optimal 
Kalman filter improves the optimal estimate of J z at the 



end of the measurement, increasing K in Eq. (38) from 
1 to 2. 



C. Single spin measurement resolution 

Single spin resolution AJ° pt < 1/2 is required for 
conditionally preparing states with spectroscopic sen- 
sitivity at the Heisenberg limit as well as the parity 
measurements needed for reading out NOON states, or 
Dicke states [57]. Single spin resolution in ensembles of 
~ 100 87 Rb atoms has recently been demonstrated using 
cavity-aided nondemolition measurements [46] . We find 
that in the far-detuned limit and with no prior knowl- 
edge, single spin resolution is reached for p < |^ , quan- 
tifying how ideal a cycling probe transition needs to be 
in order to resolve single spins. For qC ~ 1 and N ~ 10 6 , 
one would need p < 10 -7 , which is highly unrealistic for 
real multi-level level alkali-atoms due to off-resonant scat- 
teri ng fro m other hyperfine states, as will be discussed in 
Sec. VII for the case of 87 Rb. Alternatively, with high 



single atom cooperativity, qC > 8Np, single spin resolu- 
tion could be attained without a cycling transition. For 
example, with qC ~ 100 and a worst-case open transition 
with p — 1/2, single spin resolution would be reached for 
N < 25 atoms. 



VI. SPECTROSCOPIC ENHANCEMENT 

Spectroscopic sensitivity refers to the ability to resolve 
the angle through which a Bloch vector or a Dicke state 
is rotated. To first approximation, the polar angular res- 
olution is set by the conditional spin noise, discussed in 
Sec. [V] and the radius 1Z of the collective Bloch sphere 
on which the Bloch vector or Dicke state lives. In this 
section, we discuss how the radius 1Z of the collective 
Bloch sphere is reduced by the measurement due to free 
space scattering, and derive the fundamental limits to 
the spectroscopic sensitivity. The radius 1Z is propor- 
tional to the Ramsey contrast C for a CSS or a slightly 



spin-squeezed state. The radius is used here because it 
is possible to consider a conditional measurement with 
imprecision below a single spin. If all atoms remain in a 
superposition, the resulting state would be a fully sym- 
metric Dicke state or eigenstate of the operator J z , with 
1Z = N/2, but with ^ , C = 0. Nonetheless, Dicke 

states have near-Heisenberg limited spectroscopic sensi- 
tivity [58j . 

Enhanced sensitivity in one degree of angular resolu- 
tion, say the polar angle 8, can be gained at the ex- 
pense of enhanced uncertainty in an orthogonal degree 
of freedom, namely the azimuthal angle <fi. For con- 
creteness, the Bloch vector is initially prepared in a 

CSS along x with ^J^ = xN/2. The angular reso- 
lution of the polar angle for the CSS defines the SQL 
A#sql = AJ 2j css/"^ = 1/VN. If the actual angular res- 
olution is A9, then the metrolog^ically relevant squeezing 
parameter is £ m = (AOqq^/ AO) with £ m > 1 represent- 
ing a spectroscopic enhancement in sensitivity that must 
arise from entanglement. 

Angular resolution is reduced if the radius 7Z of the 
collective Bloch sphere is reduced below its initial value 
without a corresponding decrease in the spin-noise. In 
the simplest model, each free-space-scattered photon 
from an atom in a superposition state leads to the col- 
lapse of its spin into spin up or down, leading to an av- 
erage reduction in 7Z by 1/2. If the free space scattering 
rate for each spin is unchanged by the scattering process, 
then the collective Bloch sphere radius normalized to its 
initial value 7Z as a function of the number of scattered 
photons per atom m s is given by 



n 



(40) 



Note that both Rayleigh and Raman scattering lead to 
a reduction in the collective Bloch sphere radius. In cer- 
tain cases, free space Rayleigh scattering does not cre- 
ate wavefunction collapse |59| , but this requires indistin- 
guishability in the scattering process which reduces the 
information that can be extracted from the probe mode. 

Putting Eq. [36| and Eq. [40| together, the spectroscopic 
enhancement is given by 



AJl 



AJ 



z,CSS 



-2m, 



(41) 



If Raman spin-flips dominate over Rayleigh scattering 
events, the optimal spectroscopic enhancement ~ 
y qNC/p is limited by spin-flip diffusion noise as con- 
sidered in sub-section |VI A| If Rayleigh scattering dom- 
inates, then the optimal spectroscopic enhancement 

~ qNC is limited by shrinkage of the collective Bloch 
sphere radius 1Z, discussed in Sec. |VI B| The change in 
the scaling of S,^ 1 from v NC in the Raman spin-flip limit 
to NC in the cycling transition limit allows far greater 
amounts of squeezing on a cycling transition. However, 
the loss of quantum efficiency q degrades squeezing as 

~9 ona cycling transition, compared to the more 
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FIG. 17. (Color online) Optimal squeezing versus spin 
flip probability p for iVC = 1000 (black curve) and NC = 10 
(red curve) for probing in the far-detuned limit. Both curves 
take into account spin flip noise and decoherence. The curves 
assume perfect quantum efficiency q = 1 and no prior knowl- 



i.e., Cp 



co. The decoherence-limited squeezing 



£™ P = NC/e (dotted lines) is approached for p < e 2 /8NC, 
and the spin-flip-limited squeezing = */ NC/Sp (dashed 
lines) is approached for p 3> e 2 /8NC. The locus of optimal 
squeezing at the crossover point p = e 2 /8NC is 0.193/p (green 
line). 



favorable scaling £° pt 
ited regime. 



y/q in the Raman spin-flip lim- 



B. Large decoherence/cycling transition limit 



Here we consider the case where the reduction in the 
radius TZ of the collective Bloch sphere plays an impor- 
tant role in determining the angular resolution. If the 
probing is performed on a nominally closed transition 
where Raman scattering spin-flips due to probe polariza- 
tion imperfections and off-resonant scattering are very 
improbable, spin-flip diffusion noise is negligible and the 
only limit to spectroscopic enhancement is the shrinking 
of the radius TZ due to free space (Rayleigh) scattering. 
The probing is performed in the large decoherence or cy- 
cling transition regime when the optimal scattering m° pt 
that optimizes the measurement resolution of J z is not 
small. Formally, this regime occurs when the optimum 
spectroscopic enhancement calculated from Eq. ( 42 ) is 
C pt > 0.193/p. 



In the large decoherence regime, the spectroscopic en- 



hancement is given by 



-2m s 



ignoring the 



small improvement that may result from prior knowl- 
edge. An optimum is reached when the radius TZ is re- 
duced by 2 or equivalently m° pt = 1/2, yielding an opti- 
mum ££P* = l/(2eTO Pr °j). In the far-detuned limit, the 
optimal squeezing ^ rt — > qNC/e. We caution that the 
simple model presented here is not valid near the Heisen- 
berg limit because while measurement resolution below 
a single spin may be achieved, the effective spin noise 
variance must be clamped to 1/4 in order to enforce the 
Heisenberg limit on spectroscopic sensitivity. 



A. Small decoherence/spin flip limit 



VII. OPTIMAL SQUEEZING FOR 87 Rb 



Here we consider the case where the reduction in 
the radius TZ of the collective Bloch sphere may be ig- 
nored. This is justified if the optimal scattering m° pt 
that optimizes the measurement resolution of J z is small, 
m° pt <C 1, equivalently pqNC [n/T) 2 1 for probing 
in the far-detuned limit. In this regime, the radius TZ 
remains approximately 1, so that the spectroscopic en- 
hancement is primarily set by the reduction in the spin 
noise, discussed previously in Sec. W\ 



fopt 



( AJ 2 ° 
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(42) 



:-s\n 



In the 



/ AJ opt \ 2 

where ( AJ * J has been introduced in Eq. 

limit of no prior knowledge Cprcp —> oo and probing in 
the far-detuned limit, Raman spin flip noise limits the 
achievable squeezing to c^f* = y^qNC/8p. 



Using the framework above, we now analyze the lim- 
its of two separate probing schemes. From Sec. |III| we 
have shown that the measurement resolution at fixed 
free space scattering improves with cavity detuning 5 C , 
but ultimately saturates to a value set by the collective 
cooperativity parameter NC. For squeezing on a clock 
transition comprised of two hyperfine ground states, the 
maximal detuning is approximately half of the ground 
state hyperfine splitting S c = cjm/2. We are interested in 
the ultimate limits of squeezing in such a system, taking 
into account Raman spin flips and decoherence. Moti- 
vated by the fact that Raman spin flips limit the achie- 
veable squeezing on a hyperfine clock transition, we then 
analyze squeezing on a cycling transition, where Raman 
spin flips are greatly reduced, but introduce an additional 
scaling with <5 C , resulting in a region of saturation of the 
spectroscopic enhancement as NC is increased. In both 
scenarios, the ratio Whf /T of the hyperfine splitting to the 
excited state decay linewidth plays a critical role. 
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Optimal squeezing via differential measurement 
of 87 Rb clock transition 



We now consider the measurement scheme demon- 
strated by the MIT group [7] in which the pseudospin 
states were as in Ref. [T], namely the clock states of 87 Rb 
|t> = \5 2 S 1/2 ,F = 2,m F = 0) and ||) = \5 2 S 1/2 ,F = 
l,mp = 0). The bare cavity frequency is tuned to the av- 
erage of the two ground states to optically excited state 
transitions near 780 nm. This tuning ensures that an 
atom in \\.) shifts the dressed cavity resonance frequency 
by an equal, but opposite amount as an atom in |f) . The 
excited state hyperfine splitting ~ 500 MHz is much less 
than the ground hyperfine splitting and is taken to be 
zero for the following analysis. 

The problem is analyzed by extending the linearized 
two-mode model of Eq. ^ to a linearized three-mode 
model in which the atomic operator a is generalized to 
operators and d^ to yield three coupled differential 
equations, along with the same input-output relations 
given by Eq. [7] 



dt 
dt 
dt 



l -n (c) - ig (yW t (at) - v^Vj (ai)) 
(r + iw M ) (a t ) - iy/N^g (c) , 
(r - iu M ) (px) - iy/Ni~g (c) . 



(43) 



The equations are now written in a rotating frame at 
the bare cavity resonance frequency that is chosen such 
that the two optical atomic transitions are detuned by 
±Whf/2. The rate of scattering into free space is de- 
scribed by the two field amplitudes a s ,^ = VT (a^^) 
and normalized such that the rate of photons scattered 
into free space is simply M s = \a s f\ + \a s ,i\ ■ 

From the coupled set of Eq. ( 43 1 , we find that the rms 



phase shift of the transmitted light field caused by the 
rms projection noise level fluctuation in the population 
difference is: 



A<j>- 



proj 



NC 



r 



i 



i + Ncr 2 /u> 



(44) 



This expression assumes the damping rates are small 
k, T -C Whf , y/N/2 2g. The phase shift initially climbs 
with increasing atom number as y/N, but saturates to 
a maximum value A0 pr °j = \/C /2 at a critical atom 
number given by NC = (whf/T) 2 , after which the phase 
shift decreases as 1/yN. The physical interpretation 
for this decrease is that above the critical atom num- 
ber, the dressed cavity mode linewidth k' rapidly starts 
to broaden with increasing atom number. The number of 
free-space-scattered photons required to resolve the pro- 
jection noise level phase shift of the probe is 



1 



AqNC 



1 + NC 



OJhf, 



(45) 



The diffusion of the difference between the estimate of 
J z and the actual value of J z is driven by Raman tran- 
sitions that move atoms from |f) to \F = 1) or |4.) to 
\F = 2). Raman transitions between states of the same 
F (i.e., AF = 0) lead to loss of coherence, but do not 
change the coupling of the atom to the cavity mode in 
the limit where the excited state hyperfine splitting is 
neglected, as we do here. Hyperfine changing transitions 
AF ^ cause the detuning to change sign, but not mag- 
nitude, making such a process equivalent to a spin flip. 
Accounting for transition branching ratios, we find that 
to a good approximation, we can apply Eq. (36), with an 



effective spin flip probability p = 1/6. Assuming the loss 
of coherence is small, then the optimal spectroscopic en- 
hancement with respect to average probe photon number 
is 



e°pt 



y/6qNC 



1 + l.YfT 2 ,'- 



(46) 



At small N, the spectroscopic enhancement scales as 
\/6qNC, reaching a peak value of ££f* = ^3q/8uJhf /T at 
a value NC = |(oJhf/r) 2 , slightly before the maximum 
phase shift is reached. At large r NC, th e spectroscopic 
enhancement scales as = y^3q/&NC(ujhf/^) 2 - 

Taking the quantum efficiency to be q — 1, the maxi- 
mum spectroscopic enhancement for 87 Rb is quite large 
at 28 dB. The exact details of the full measurement 
sequence (i.e., whether rotations such as 7r-pulses are 
used to cancel sources of technical noise as was done in 
Refs. Q][2]) are needed to construct an optimal estima- 
tor of J z , but, at best, a 3 dB further improvement may 
result. 

Because C does not depend on the cavity length, the 
optimum N for peak spectroscopic enhancement scales as 
(u>Q/-F)(u;hf/r) 2 , where wq is the cavity mode waist, and 
F is the cavity finesse. More fundamentally, no change in 
the cavity geometry (wo and I) or finesse F changes the 
maximum obtainable enhancement in spectroscopic sen- 
sitivity. This enhancement is determined solely by the 
atomic properties. Figure. [18] shows the spectroscopic 
enhancement versus atom number for a range of techno- 
logically feasible cavity finesses. 

Resolving very small phase deviations or small fre- 
quency shifts imposes technical challenges that are mod- 
ified with cavity geometry or finesse, as shown by the 



probe frequency (Fig. [19| and probe phase shift (Fig. 20 ) 



resolutions required to obtain the spectroscopic sensitivi- 
ties shown in Fig. [18] All three figures assume the cavity 
geometry of Ref. pQ, I = 1.91 cm and wq = 71 fim. 



Finally, we note that Fig. 19 shows the technical require- 
ment on probe frequency resolution is more relaxed above 
the optimum N compared to achieving the same spectro- 
scopic enhancement at a value below the optimum N. 
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FIG. 18. (Color online) The fundamental optimum spectro- 
scopic enhancement for differential probing of the 87 Rb clock 
transition, as performed in Ref. [7]. The calculations assume 
a net quantum efficiency of q = 1 and the cavity geometry 
described in the text and Ref. pQ. The purple, light blue, 
green, orange, and red curves correspond to cavity finesses of 
F — 10 2 , 10 3 , 10 4 , 10 5 , and 10 6 , respectively. All these finesses 
are experimentally feasible. The atom number N at which 
the spectroscopic enhancement is maximized scales with the 
cavity mode waist wo and cavity finesse F as Wq/F. 



FIG. 20. (Color online) Phase resolution with which the 
transmitted probe light must be measured to obtain the 
spectroscopic enhancements shown in (and under the same 
conditions as) 




10 10 10 10 10 10 10 10 10 10 10 
Atom Number N 

FIG. 19. (Color online) Frequency resolution with which the 
relative frequency of the probe and dressed cavity mode must 
be measured to obtain the spectroscopic enh anc ements shown 
in (and under the same conditions as) Fig. 18 . The purple, 



light blue, green, orange, and red curves correspond to cavity 
finesses of F = 10 2 , 10 3 , 10 4 , 10 5 , and 10 6 , respectively. 



B. Optimal squeezing on Rb cycling transition 

Having seen that the optimal squeezing on a clock tran- 
sition is fundamentally limited by Raman spin flips, we 
consider a situation in which Raman spin flips are re- 
duced, namely probing on a cycling transition, and show 
that larger amounts of squeezing are possible in this con- 
figuration than on the clock transition [46 [ 160] . 

As a concrete example of how a cycling transition 
can be used to enhance probing, we consider the cy- 
cling transition in 87 Rb, \-\) = \F = 2,m F = 2) to 
|e) = \F = 3', rap = 3) at wavelength 780 nm. The spin 
down state is chosen as \l) = \F = l.mp = 1). For the 



Fig. 18 



The purple, light blue, green, or- 
ange, and red curves correspond to cavity finesses of F = 
" , 10 3 , 10 4 , 10 s , and 10 6 , respectively. 



following, the probing scheme with relevant energy levels, 
dipole matrix elements, decay branching ratios, dressed 
mode frequencies, and probe laser detunings are shown 
and defined in Fig. |21| Here we will extend the previous 
models of the precision of the estimation of J z and the 
loss of signal due to wavefunction collapse to capture the 
essential physics for this system. Key results are that 
there exists a region of saturation, or universal spectro- 
scopic enhancement, set only by atomic properties and in 
which varying atom number and cavity finesse can have 
little impact. However, unlike in the previous section, the 
asymmetry in the cavity coupling to |f) and \\) allows 
this saturation region to be surpassed at large values of 
NC. 

We must first consider what limits the rate of Raman 
scattering processes that can lead to diffusion of the spin 
projection J z . The probe polarization can be set to pure 
a + to better than 10 -4 , so that Raman scattering from 
|t) is suppressed to at least this level or greater. The more 
fundamental Raman scattering limitation arises from the 
finite hyperfine splitting cjhf = 27r x (6834 MHz). Specif- 
ically, atoms in | J.) can non-resonantly Raman scatter 
probe photons from |e') = \F = 2', nip = 2). 

In the following discussion, the quantity m s is impor- 
tantly defined as the average number of probe photons 
(normalized to the total atom number) Rayleigh scat- 
tered into free space by atoms in 1 1) • All other scat- 
tering processes will be scaled from this quantity using 
the quantities defined in Fig. |21| The key parameters for 
rescaling are the ratio of the dipole matrix elements r — 
Mle'/M^e — l-v/2, the decay branching ratio B 3e > = 1/6 
from |e') to |3) = \F = 2,m F = 1), and the detunings of 
the probe light S e and 8 e > = 5 e — Whf — w ohf from resonance 
with the transitions |f) — > |e) and ||) — > |e') respectively. 
As in the previous section, we neglect the excited state 
hyperfine splittings so that S e i « 5 e — cjhf . 

The rms imprecision in the estimate of J z relative to 
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FIG. 21. (Color online) Probing scheme for the cycling tran- 
sition in 87 Rb. The ground hyperfine states F — 2 and F = 1 
are split by ojhf = 27r x (6834 MHz). The values of m_F are la- 
beled across the top. The relevant F' = 3 and F' = 2 excited 
D2 transitions at wavelength 780 nm are also shown with the 
excited state splitting w e hf = 2n x (267 MHz). The pseudo 
spin-1/2 system is here composed of |t) = \F = 2,mp = 2) 
and |4-) = \F = l,raj? = 1). Ideally, the cr + -polarized probing 
laser couples only |f) to the optically excited state |e) with 
dipole matrix element Mf e at a frequency detuning 8 e that 
is approximately equal to the dressed-cavity mode frequency 
uj-. In this illustration, S e is negative. By dipole selection 
rules, \e) can only decay back to |f). However, the same 
probe laser also couples | \) to the single excited state |e') 
with dipole matrix element Mj, e / and larger detuning from 
resonance S e r = 5? — uihf — oj c hf • The ratio of matrix elements 
is M± e i/Mf e = 1/V2. Finally, the state |e') decays to states 
|4-), |t), and 1 3) with fractional branching ratios B± e i = 1/2, 



B 



1/3, and B 3e i =1/6 respectively. 



the projection noise level can be approximately modeled 
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Ap^m s + p 3 m s 



(47) 



Starting in order of physical significance, the first and 
second terms arise from diffusion of J z caused by Raman 
transitions from \l) — > |t) and |3) with effective proba- 
bilities p^ and j»3 given approximately by 



B tl3 cr 2 



5,, 



(48) 



In this simple treatment, Raman decays to |3) are treated 
as loss, as reflected in the smaller numerical pre-factor in 
front of the second term of Eq. ( 47 ) . 



The third term in Eq. (47) is modified to reflect that 



both states can interact with the probe at large detun- 
ings such that the dressed cavity mode frequency is less 
sensitive to quantum projection noise in J 2 , and thus 
more probe photons must be used to resolve J z at the 



projection noise level, i.e. 



R 



(49) 



ray 



Here, is defined by Eqs. (Q, and @. Indis- 

tinguishability is accounted for by 



i?ray=(l-r- 



(50) 



1 



Note that i? ray < 1 with an asymptotic value of -R ray 
at large detunings. 

There are two effects that are neglected in Eq. [49] by 
first assuming they are small, and then verifying this to 
be the case after the calculations. First, in applying the 
dressed cavity linewidth result for k 1 from Eq. ([8]), we 
assume that the cavity mode is negligibly further broad- 
ened by atoms in state | \) . By estimating the addi- 
tional broadening evaluated at the optimal cavity detun- 
ing and average number of scattered photons, we find 
that the optimal spectroscopic enhancements calculated 
in Fig. [22] are reduced by < 0.3 dB due to the neglected 
mode broadening. Second, we assume that the dressed 
mode frequency w_ calculated from Eq. Q is only modi- 
fied by a small fraction by atoms in state \\.) . Again, this 
assumption is verified to be the case at the optimal cavity 
detuning and the average number of scattered photons, 
with the exception of the case where N > 10 s and cavity 
finesse F = 100, as shown in Fig. [22] where several dB of 
deviations are possible due to this effect. 

Next, we consider how collapse due to free space scat- 
tering reduces the radius of the collective Bloch sphere, 
specifically 



n 



(51) 



where the partial cancellation of wavefunction collapse 
due to indistinguishable Rayleigh scattering off of both 
|t) and ||) (see [59]) is accounted for by J2 ray . 

The term i? ram accounts for Raman scattering from |4_) 
to |3). 



Rr. 



B 



3e" 



(52) 



As before, we assume that Raman scattering to state |3) 
is equivalent to atom loss. Note also that i? ra m < B^r 2 . 

Equations (47)-(52) are used to numerically estimate 
the optimal spectroscopic enhancement shown versus 
atom number N in Fig. [22] for a range of technologically 
reasonable cavity finesses and assuming the cavity geom- 
etry of Ref. pP (cavity length I = 1.91 cm, mode waist 
size wq = 71 /xm), and perfect quantum efficiency q = 1. 
The optimization is done with respect to both m s and 
the dressed cavity mode frequency w_ (tuned by chang- 
ing the bare cavity frequency). The mode frequency u;_ 
at the optimum is shown in Fig. [23] The loss of signal 
due to wavefunction collapse and scattering to |3) at the 



optimum is shown in Fig. 24 
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FIG. 22. (Color online) Theoretical optimal spectroscopic en- 
hancement f££ (solid curves) in dB relative to the SQL versus 
total 87 Rb atom number N for the current cavity geometry of 
Ref. [T] and probing near the D2 cycling transition as shown in 
Fig. |21| The optimization is performed for a range of techni- 
cally realizable cavity finesses F = 10 2 , 10 3 , 10 4 , 10 5 , and 10 6 
(blue, cyan, green, yellow, and red curves respectively) under 
the assumption of perfect quantum efficiency q = 1. Note 
the unphysical regime below the Heisenberg limit (disallowed 
grey region). A region of saturation of £J)f versus N occurs 
near ss ojhf/r shown by the labeled horizontal dashed 
line. The physical origin for the saturation region arises from 
competition between the scaling of the off-resonance Raman 
scattering probabilities and the dressed cavity mode broad- 
ening, as described in the text. The points of 3 dB deviation 
from NC scaling at low atom number (solid circles) and y/NC 
scaling at large atom number (open circles) occur when the 
solid curves cross ££f* = 27 dB and 33 dB, respectively. The 
approximate transition from q to ^fq scaling occurs when the 
solid curves cross ^ l = 28 dB, the point at which the scaling 
is approximately q 3 ^ 4 . 
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FIG. 23. (Color online) The dressed mode frequency ui- that 
optimizes the spectroscopic enhancement £ m for the condi- 
tions described in Fig. |22[ where again the cavity finesses F = 
10 2 , 10 3 , 10 4 , 10 s , and 10 6 correspond to blue, cyan, green, 
yellow, and red curves, respectively. The open and closed 
circles indicate the locations of the 3 dB points in Fig. [22] 
The mode frequency is normalized to the hyperfine splitting 
w M /27r = 6.834 GHz. Note that uj- « F at LU-/ui hi « 10~ 3 . 
Comparing to Fig. |22| one finds that the transition to spec- 
troscopic enhancement scaling as ££f oc yfqNC occurs near 
u)- ~ 1.8 x cjhf. Above this point, the spin flip probabilities 
change little with w_ . 



FIG. 24. (Color online) The scattering m 3 that optimizes the 
spectroscopic enhancement ££f * shown in Fig. |22| Again, the 
cavity finesses F = 10 2 , 10 3 , 10 4 , 10 5 , and 10 6 correspond to 
the solid blue, cyan, green, yellow, and red curves, respec- 
tively. The open and closed circles indicate the locations of 
the 3 dB points in Fig. |22| Comparing to Fig. |22| one finds 
that the transition to spectroscopic enhancement scaling as 
££f oc qNC occurs when m s ~ 0.25, indicating that at larger 
N, far off-resonance Raman scattering begins to create signif- 
icant diffusion of J z that limits the precision of the estimation 
of J, as discussed in the main text. 



At low atom number, the spectroscopic enhancement 
scales as ~ qNC. At high atom number, the spec- 
troscopic enhancement scales as ~ y/qNC. There 
is an intermediate region of atom number for which the 
spectroscopic enhancement is relatively flat versus atom 
number with ££P' ~ (jy^/T. The physical origin of this 
plateau arises from the form of the critical detuning of 
Eq. pi. 



In the good cavity limit, the scattering necessary to 
reach projection noise level sensitivity falls as l/(5 4 

for S c < 5°, making it beneficial to operate with |cj_| ~ 
S c > r 'a/ ' NC I 'q. However, the Raman transition proba- 
bilities continue to grow quadratically with the de- 
tuning, while detuning farther no longer rapidly reduces 

Assuming the critical detuning 5° is optimal for the 
reasons above, and in the limit of < uj^f, the Raman 
transition probabilities scale as p^, 3 ~ (NC '/ 'q)(T /whf) 2 , 
while mP r °j ~ 1/qNC. Optimizing the total noise in 
our estimate of J z using Eq. (47) with respect to m s 
reproduces the observed plateau value c^j 3 * ~ u>hf/T ~ 
10 3 . The plateau region is exited at low atom number 



when loss of signal [described by Eq. (51 )] dominates the 
reduction in spectroscopic enhancement, as illustrated by 
the loss of signal due to wavefunction collapse shown in 
Fig. 
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At high atom number, the plateau region is 
exited when the optimum mode frequency becomes large 
compared to the hyperfine splitting |a;_| > Wht, as shown 
in Fig. [24} 

Importantly, this analysis shows that there is a range 
in which increasing either finesse or atom number can 
have little effect on the optimal spectroscopic enhance- 
ment achieved. Also, note that the value of the plateau 
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does not depend on the cavity geometry, and therefore 
represents a universal value that depends only on the 
atomic properties. One can also show that the range 
of atom number spanned by the plateau and the posi- 
tion of the first corner scales roughly as N ~ Wq/F, but 
does not depend on the cavity length I. Finally, for atom 
numbers below 10 3 , it appears possible to both prepare 
and readout states near the Heisenberg limit using this 
approach and technologically feasible cavity finesses. In- 
deed, Ref. [46] recently demonstrated single atom mea- 
surement resolution for N ~ 100 using the approach de- 
scribed here. 



VIII. COHERENT REINITIALIZATION FOR 
REDUCED RECOIL HEATING 

In this final brief section, we consider whether coher- 
ence preserving quantum non-demolition measurements 
can also be used to reduce the total photon-recoil heat- 
ing in a Ramsey spectroscopy measurement cycle. A low 
number of photon recoils imparted in a complete exper- 
imental cycle, not just the readout, would allow neutral 
atom traps to operate more like ion traps, with their 
advantages of high duty cycle and high repetition rates. 
We use the canonical example of Ramsey spectroscopy to 
show how low photon recoil in a complete experimental 
cycle can be achieved. Imagine reading out the Bloch vec- 
tor projection J z at the end of a Ramsey sequence with a 
QND measurement at or below the projection noise limit. 
If one were to repeat the experimental cycle by optically 
repumping all the atoms back into each atom in the 
ensemble would experience ~ s/2 photons worth of re- 
coil heating, where s is the average number of photons 
scattered by an atom to get from |f) to \l) during the 
repumping process. In Fig. |25| we show how recoil heat- 
ing can be reduced by coherently reinitializing most of 
the atoms into ||.) using a microwave 7r/2 pulse before 
turning on the repumping. 

The loss of contrast e due to dephasing during the 
Ramsey evolution time is assumed to be small for the 
results presented in this section to be valid. Dephasing 
causes the off-diagonal elements of the density matrix 
representing the state of the ensemble to decay, resulting 
in a mixed state. The mixed state can be regarded as a 
tensor product of a pure state representing a smaller en- 
semble of size N(l — e) atoms in a coherent superposition 
of |t) and ||) [Fig. |25[a)], with pure states representing 
the decohered atoms with eN/2 atoms collapsed in |f) 
and | J.) on average. The QND readout causes an ad- 
ditional (1 — p)M s atoms to collapse into |f) and pM s 
atoms to collapse into ||), provided m s — M s /N <C 1. 
Following a microwave ir/2 pulse that rotates most of the 
atoms into |4-), we find that approximately s(m s + e)/2 
photon recoils per atom is required to reinitialize the en- 
tire ensemble into The contribution from quantum 
backaction into J y is of order sm s /NmP r ° 3 and may be 
safely ignored for small spectroscopic enhancements. 



(a) Atomic spin state after QND measurement: 



A.I 



y, QBA 



N — M s — eN atoms 



{l-p)M s + — pM s + — 



atoms 



atoms 



(b) After 7i/2-rotation about y: 




M s +eN 
(l\f) = atoms 



(c) After repumping: 




FIG. 25. (Color online) Coherent reinitialization scheme to 
reduce photon recoil heating in a complete experimental cycle, 
(a) Dephasing during the Ramsey evolution can be described 
as a shrinking of the coherent part of the Bloch vector from 
length N/2 to length N(l — e), where e<1 characterizes the 
amount of dephasing. The dephased atoms can be represented 
as if they had collapsed into |t) and |4-) with eN/2 atoms in 
each state. Performing a QND readout with m s = M a /N <C 1 
results in further shrinkage of the coherent Bloch vector and 
additional collapse of the atoms, (b) Performing a microwave 
7r/2 pulse about y rotates the coherent Bloch vector to the 
south pole and the collapsed atoms onto the equator of the 
Bloch sphere, (c) Applying optical repumping beams reini- 
tializes the ensemble to the south pole, with the photon recoil 
heating arising from repumping the atoms in superposition 
states into |4-). The contribution arising from the quantum 
backaction of repumping back into a CSS is negligible pro- 
vided the spectroscopic enhancement is small compared to 
the Heisenberg limit. 



In summary, we have shown that it is possible to im- 
part only s(m s + e)/2 optical photon recoils per atom 
through the coherent reinitialization scheme, as opposed 
to the much larger s/2 recoils per atom without coherent 
reinitialization. Note the reason why the scheme works 
is because we have traded optical photon recoils with mi- 
crowave photon recoils, which impart recoil energy ~ 8 
orders of magnitude smaller than that of optical photons. 
If an optical transition separated [f) from ||,), such as in 
optical lattice clocks, the scheme presented would not 
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reduce recoil heating substantially. 

IX. CONCLUSIONS 

In conclusion, we have presented detailed expressions 
for how cavity-aided nondemolition measurements of 
atomic populations scale with key experimental param- 
eters: cavity linewidth, cavity geometry, collective coop- 
erativity, and Raman transition probabilties. This anal- 
ysis will be important in future experiments for reduc- 
ing photon-recoil heating and for producing highly en- 
tangled states with enhanced spectroscopic sensitivities 
relative to the SQL. We have given experimental details 
for two different probing schemes in 87 Rb and estimated 
fundamental limits on conditional spin-squeezing in en- 



sembles of Rb atoms. We expect that this analysis will 
guide future experimental efforts on squeezing on a cy- 
cling transition, specifically in 87 Rb and can be straight- 
forwardly adapted to other atoms with ground state hy- 
pcrfine structure such as Cs and other alkali elements 
used in spectroscopy and quantum information experi- 
ments. 
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